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s  Umax 


An  electron  bean  probe  was  used  to  obtain  measurements  of  rota¬ 
tional  temperature  and  gas  density  in  supersonic  nitrogen  jets  expanding  froa 
room  temperature .  The  parameter  PQd  (where  P0  is  the  stagnation  pressure  in 
torr,  and  d  is  the  orifice  exit  diameter  in  xma)  was  varied  froa  15  to  UQo 
torr-nm.  This  corresponds  to  a  sonic  Reynolds  cumber  range  of  230  to  9,320. 
Density  measurements  were  made  using  a  photomultiplier  with  an  interference 
filter  centered  at  3300  A°  and  having  a  half -width  of  100  A°-  The  experimental 
density  data  follow  the  axial  isentropic  density  distribution  in  regions  of 
rotational  nonequilibriua.  Eotatioual  temperature  measurements  were  obtained 
froa  rotational  spectra  of  the  first  negative  0-0  band  at  391^  A°  using  a 
high  dispersion  spectr  jgraph.  The  experimentally  detemiced  rotational  tem¬ 
perature  values  folio*  the  axial  isentropic  temperature  distribution  in  a  free 
jet  down  to  auproxim?  „ely  85  °K  for  PQd  =  15  torr-ra,  and  50  °K  for  PQd  =  ^80 
torr-mm,  with  a  scatter  of  about  t  3p-  Below  these  temperatures,  the  Tg  data 
depart  froa  the  isentropic  curve  and  freeze  at  a  constant  temperature,  which 
is  dependent  on  the  value  cf  ?Qd. 

A  shock  holder  was  inserted  in  the  jet  and  a  number  of  shock 
waves  in  the  range  M  =  h  to  M  =  15  were  investigated.  Density  profiles  through 
the  shock  waves  were  obtained.  Botational  spectra  indicate  a  large  departure 
from  a  3oltzmann  distribution  in  the  rotational  levels  in  the  center  of  a  shock 
front.  This  effect  is  small  at  U  -  h,  but  very  pronounced  at  M  =  15.  An 
apparent  non-Boltzmann  rotational  distribution  in  the  jet  expansion  flow  was 
also  observed. 
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I.  IKTROBUCTOBY  DISCUSSION  „  - 

The  unconfined  expansion  from.  a  £onic  orifice  into  a  low  pressure 
chaster  (referred  to  as  an  underexpanded  or  free  jet)  has  proven  to  be  -a  very 
popular  and  useful  source  of  a  high  velocity,  low  density  gas  flow. 

As  outlined  in  the  recent  review  paper  by  French  (Ref.  l),  this 
basic  expansion  flow  is  especially  useful  for  facilities  having  a  limited  pump¬ 
ing  capacity,  since  the  over exp ansi on  in  the  jet  brings  the  flow  to  exceedingly 
high  local  Each  numbers.  The  availability  of  theoretical  treatments  based  oh 
isentropic  assumptions  enables  flow  properties  to  be  computed  which  are  applies 
able  over  a  wide  range  of  conditions  of  experimental  interest. 

A  number  of  authors  in  recent  years  have  investigated  ufceoretieally 
the  flow  field  of  an  underexpanded  jet.  Oven  and  Thornhill  (Ref.  ?)  if  ere  the 
first  to  consider  the  problem  and  used  a  method  of  characteristics  solution  for 
a  gas  expanding  into  a  vacuum.  When  the  jet  expands  into  a  region  of  finite 
pressure,  however,  a  complicated  flow  pattern,  such  as  diagrammed  in  Fig,  1 
develops.  A  method  of  characteristics  solution  will  be  valid  up  to  the  first 
Mach  disc,  since  the  flow  within  the  barrel  shock  is  uninfluenced  by  pressure 
changes  along  the  jet  boundary.  Investigations  at  NASA  by  love  and  Associates 
(Bef.  3)  and  Vick  et  al  (Bef.  4)  have  dealt  with  both  theoretical  and  experi¬ 
mental  studies  of  free  jets,  where  the  position  of  the  barrel  shock,  jet 
boundary  and  Mach  disc  was  computed  for  a  variety  of  conditions. 

Extended  method-of-eharacteristics  calculations  by  Woolf  for 
high  stagnation  to  ambient  pressure  ratios  are  presented  in  Bef.  5-  Tbs  center- 
line  Mach  number  distributions  obtained  from  these  calculations  are  shown  in 
Fig.  2  for  two  values  of  specific  heat  ratios  ( 7 ) .  Recent  theoretical  and 
experimental  free  jet  studies  are  summarized  in  the  excellent  reference  paper 
by  Ashkenas  and  Sherman  (Ref.  6),  who  have  shown  good  agreement  between  experi¬ 
ment  and  theory. 

Relaxation  processes  in  expanding  gas  flows,  particularly  flows 
occurring  in  rocket  and  shock  tunnel  nozzles,  have  been  studied  both  theoretically 
and  experimentally  for  some  time.  Chemically  reacting  expansions,  for  example, 
have  bean  investigated  by  Hall,  et  al  (Bef.  7),  while  the  problem  of  relaxation 
of  the  vibrational  degrees  of  freedom  of  diatomic  molecules  has  been  studied 
theoretically  by  Stollery  and  Smith  (Ref.  8),  and  Bray  (Ref.  9)*  The  calcula¬ 
tions  of  Treanor  (Ref.  10)  have  shown  the  possibility  of  a  departure  from  a 
Boltzmann  population  distribution  in  the  vibrational  levels  of  a  diatomic 
molecule  during  a  rapid  expansion.  Recent  experimental  investigations  such  as 
those  of  Hurle,  Russo  and  Hall  (Ref.  li),  and  Nagaaatsu  and  Sheer  (Ref.  12), 
for  example,  have  developed  techniques  to  study  relaxing  expanding  gas  flows. 

All  of  the  above-mentioned  studies  have  been  at  elevated  stagna¬ 
tion  temperatures,  greater  than  2000°K,  where  there  is  sufficient  excitation  of 
the  internal  vibrational  degrees  of  freedom.  The  results  of  all  indicate  a 
departure  from  equilibrium  as  a  vibrating  diatomic  gas  undergoes  a  sufficiently 
rapid  expansion,  hue  to  the  finite  reaction  rates  involved  (i.e.  the  rather 
large  number  of  collisions  needed,  1000-10,000  to  maintain  equilibrium  among 
the  vibrational  and  translational  degrees  of  freedom),  the  slower  vibrational 
modes  are  unable  to  follow  translation,  and  depart  from  equilibrium  in  the  ex¬ 
pansion.  If  the  population  in  the  vibrational  levels  remains  Boltzmann,  then 
the  results  show  that  the  vibrational  temperature,  Ty,  peels  away  from  the 
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translational  temperature,  T^.  Finally,  sufficiently  far  downstream,  collisions 
become  so  rare  that  the  vibrational  temperature  change  is  very  small,  and.  it  is 
said  to  be  nearly  frozen  at  a  value  higher  than  the  local  translational  tempera¬ 
ture. 

Since  most  free  jets  under  study  are  at  room  temperature  stagn a» 
tion  conditions,  there  is  no  vibrational  excitation  of  importance  in  the  case  of 
nitrogen.  Also,  owing  to  the  extreme  rates  of  expansion,  the  vibrational  modes 
are  frozen  immediately  at  their  room  temperature  (ground  level)  population. 
However ,  in  tht  expansion  extremely  low  temperatures  are  reached,  as  can  be 
seen  from  Fig.  3?  where  the  temperature  distribution  is  plotted  corresponding 
to  the  Mach  number  distribution  of  Fig.  2,  and  at  sufficiently  low  density  one 
may  expect  the  rotational  degree  of  freedom  to  go  out  of  equilibrium  with  the 
translational  degrees  of  freedom.  Although  the  number  of  collisions  needed 
to  equilibrate  rotation  with  translation  are  small  (approx.  5  to  10),  the  low 
densities  encountered  in  such  flows  may  lead  to  conditions  where  the  number  of 
collisions  are  not  sufficient  to  maintain  even  this  fast  internal  degree  of 
freedom  in  equilibrium  with  translation.  A  recent  theoretical  approach  to  this 
problem  is  that  of  Knuth  (Ref.  13),  in  which  the  axial  distance  from  the  orifice 
where  rotational  freezing  may  occur  is  predicted,  based  on  a  room  temperature 
(300°K)  collision  number. 

In  this  type  of  calculation,  the  rotational  temperature  is 
assumed  zo  be  in  equilibrium  with  the  translational  temperature  until  a 
sufficiently  low  temperature  and  density  (a  function  of  stagnation  conditions) 
are  reached  in  the  expanding  jet.  At  this  point,  the  rotational  temperature 
is  assumed  to  depart  immediately  and  remain  frozen  at  a  constant  value,  while  the 
translational  temperature  continues  to  decrease  in  the  expansion.  This  approach 
is,  at  best,  only  ajjproximate  for  two  reasons:  the  rotational  collision  number 
at  the  low  temperatures  (~  20°K)  encountered  in  free  jets  is  unknown,  and 
secondly,  the  rotational  degree  of  freedom  dees  not  suddenly  freeze  at  a  con¬ 
stant  temperature  T^,  but,  may  depart  gradually  from  the  translational  degree  of 
freedom.  However,  an  analysis  ol  this  type  given  by  Knuth  (Ref.  13)  does  serve 
to  delineate  operating  regimes  of  stagnation  pressure  (PG)  and  orifice  exit 
diameter  (a)  under  which  rotational  nonequilibrium  may  occur  in  an  underexpanded 
jet. 

Experimental  determination  of  rotational  nonequilibrium  in  an 
underexpanded  jet  has  proven  to  be  quite  elusive.  Impact  probe  measurements 
such  as  those  reported  by  Ashkenas  and  Sherman  (Ref.  6),  for  example,  will  not 
define  this  departure  from  equilibrium  for  most  jets  under  study.  The  impact 
probe  responds  essentially  to  the  quantity  (p  u2),  where  p  is  the  density  and 
u  is  the  velocity  in  the  jet.  The  flow  velocity  in  an  underexppded  jet 
approaches  the  limiting  velocity  in  an  expansion,  (y  +  l/y  -  1)2  a*,  (where  a* 
is  the  sonic  velocity  at  M  =  l) .  The  density  then,  falls  off  approximately  as 
(x/d)-^,  sinie  the  jet  can  be  approximated  as  a  radial  expansion  from  a  source 
near  the  orifice  (Ref.  6).  In  addition,  for  the  usual  experimental  conditions, 
the  rotational  energy  may  be  excited  in  the  shock  wave  ahead  of  the  impact 
probe,  further  adding  to  the  complications.  As  pointed  out  by  Knuth  (Ref.  13), 
the  difference  in  inpact  pressure  at  a  given  x/d  is  only  about  20%  for  the  cases 
cf  immediate  rotational  freezing  upstream  of  the  orifice,  and  the  freezing  of 
rotational  degrees  of  freedom  at  a  high  Mach  number  downstream  of  the  orifice. 

There  is,  however,  a  technique  'which  may  be  used  to  investigate 
the  rotational  temperature  distribution  in  an  underexpanded  jet,  to  determine 


the  departure,  if  any,  of  the  rotational  degree  of  freedom  from  equilibrium. 

This  is  the  electron  beam  fluorescence  technique  described  originally  by  Muntz 
(Ref.  1*0-  lbs  usefulness  has  been  further  demonstrated  by  Muntz  and  Marsden 
(Ref.  15),  Marsden  (Ref.  16),  Muntz,  Abel  and  Maguire  (Ref.  17),  and  Sebacher 
(Ref.  18) .  Recently,  this  technique  has  been  used  by  a  number  of  investigators 
to  measure  rotational  temperatures  in  expanding  flows.  Sebacher  and  Duckett 
(Ref.  19)  and  Petrie  (Ref.  20),  for  example,  have  obtained  measurements  in 
high  temperature  expansions  with  an  arc-tunnel  nozzle.  In  addition,  Rbbben 
and  Talbot  (Ref.  21)  have  investigated  underexpanded  nitrogen  sets  expanding 
from  room  temperature.  _  - 

The  predominant  contribution  to  -the  fluorescence  observed  in 
nitrogen  at  low  densities  arises  from  the  first  negative  system  which  repre¬ 
sents  radiation  from  excited  molecular  nitrogen  ions.  Muntz  presented  a  theo¬ 
retical  analysis  of  the  rotational  intensities  in  the  bands  of  this  system, 
and  found  experimentally  that  the  rotational  temperature,  Tp,  (and  in  more 
detail  the  relative  population  distribution  in  the  various  rotational  energy 
levels)  of  the  nitrogen  molecule  before  electron  excitation  could  be  obtained 
by  using  this  theory.  It  is  this  electron  beam  fluorescence  technique  that 
has  been  used  in  the  present  investigation  to  study  departures  from  rotational 
equilibrium  in  an  underexpanded  jet. 

Since  a  free  jet  represents  a  convenient  source  of  high  Mach 
number  flow,  it  may  be  used  in  the  study  of  the  structure  of  strong  shock  waves 
by  inserting  a  shock  holder  in  the  high  velocity  flow  ahead  of  the  first  Mach 
disc.  The  electron  beam  technique  in  this  system  will  yield  data  on  the 
rotational  distribution  within  a  shock  wave.  In  addition,  by  measuring  the 
intensity  of  the  total  fluorescence  due  to  the  beam,  the  density  distribution  in 
a  shock  wave,  and  thus  a  measure  of  the  shock  thickness,  can  be  obtained. 

Earlier  studies  of  shock  waves  have  utilized  low  density  wind  tunnels,  such  as 
that  used  by  Sherman  (Ref.  22)  in  his  heated  wire  experiment,  and  shock  tubes. 
Shock  thickness  measurements  have  been  obtained  by  a  variety  of  techniques  in 
shock  tubes,  for  example,  the  optical  refractivity  method  employed  at  Princeton 
by  Linzer  and  Hornig  (Ref.  23).  The  attenuation  of  an  electron  beam  was  used 
by  Ballard  and  Venable  (Ref.  24),  Russell  (Ref.  25),  Schultz-Grunow  and  Prohn 
(Ref.  26),  while  the  eomprehens ive  measurements  of  shock  thickness  of  Camac 
(Ref.  27)  utilized  the  electrons  scattered  from  the  primary  beam.  Camac' s 
technique  yields  directly  a  density  trace  through  the  shock  wave.  All  of  the 
above  mentioned  methods  depend  on  the  ratio  of  the  density  behind  the  shock 
to  the  density  ahead  of  the  shock  to  determine  if  rotational  equilibration  has 
taken  place  within  the  shock  front.  However,  by  observing  the  rotational  spectra 
arising  from  beam  fluorescence,  it  can  be  determined  whether  or  not  there  is 
a  departure  from  rotational  equilibrium  within  the  shock  front  itself.  It  may 
be  expected,  for  example,  that  non-Boltzmann  rotational  populations  may  occur 
in  the  centre  of  shock  waves.  The  electron  beam  technique  will  indicate 
whether  or -not  this  is  in  fact  so. 

The  present  investigation,  utilized  the  electron  beam  fluorescence 
to  study  the  state  of  the  rotational  degrees  of  freedom  in  underexpanded  jets 
and  shock  waves.  Section  II  of  this  report  discusses  the  theory  of  rotational 
temperature  measurements  with  the  electron  beam  and  the  experimental  apparatus 
is  described  in  Section  III.  The  final  results  and  discussions  are  presented 
in  Sections  17  and  V. 


II.  ROTATIONAL  TSMESRATDRS  MEASUREMENTS  USING  AN  ELECTRON  BEAM 


A  beam  of  energetic  electrons  when  passing  through  .a  low  density 
gas  excites  a  fluorescence  that  can  be  used  to  obtain  information  concerning  - 
the  gas  molecules  prior  to  their  excitation.  The  two  properties  of  interest 
in  the  present  investigation  are  the  rotational  temperature.  Is,  of  the  nitro¬ 
gen  molecule  before  excitation,  and  the  density  of  the  gas.  If  the  electron 
beam  is  well  collimated,  the  fluorescence  is  maintained  in  a  cylindrical  region 
around  the  heam.  With  a  carefully  designed  optical  system,  it  is  possible  to 
select  a  small  elemental  volume  or  "point"  anywhere  along  the  beam  length  to 
obtain  emission  data.  As  pointed  out  by  Muntz,  et  al  (Ref.  17),  for  this 
technique  to  be  useful,  at  least  one  .relatively  strong  emission  line,  must  be 
observed  when  the  gas  is  bombarded  by  energetic  electrons.  Bor  nitrogen 
at  low  pressure  (i.e.  less^than  1  mm  Rg.),  the  predominant  radiation  is.  emitted 
from  transitions  of  the  jSg'  molecular  ion  (the  1st  negative  system)  with  the 
0-0  vibrational  transition  being  tbe  strongest  emitter  (see  Ref .  15).  As  the 
pressure  is  increased,  the  main  contribution  to  the  fluorescence  arises  from 
the  N2  second  positive  system.  Ibis  effect  was  exhibited  in  the  results  of 
Davidson  and  0‘Neil  (Ref.  28),  where  the  spectral  intensity  of  the  beam  fluore¬ 
scence  was  observed  at  pressures  hanging  from  5  him  2g.  up  to  approximately  one 
atmosphere. 


In  addition,  the  excited  state  lifetime  should  he  less  than  10“7 
seconds,  to  eliminate  the  possibility  of  the  fluorescence  being  carried  down¬ 
stream  in  a  high  'velocity  gas  flow.  Bor  the  nitrogen  first  negate  :e  system, 
and  in  particular  the  0-0  transition,  the  lifetime  is  approximately  8  x  10“° 
seconds,  Nicholls  (Ref.  29).  Thus,  at  a  flow  velocity  of  20,000  feet  per  second, 
an  emitting  particle  would  he  carried  downstream  approximately  0.4  mm.  Bor  .a 
heam  width  of  one  mm,  this  drift  would  male  good  resolution  of  the  beam  center- 
line  very  difficult,  however,  the  limiting  velocity  obtained  in  a  free  jet 
expansion  from  room  temperature  is  approximately  2,500  feet  per  second  for  a 
7  =  1.4  gas.  The  drift  in  this  case  is  .about  0.05  him,  and  it  is  possible  to 
obtain  better  spatial  resolution  by  observing  the  center  portion  of  the  beam. 

Since  the  beam  fluorescence  arises  from  a  direct  excitation  - 
spontaneous  emission  sequence,  there  is  an  upper  bound  in  gas  density  which 
indicates  the  onset  of  ccllisional  quenching  (i.e.  tbe  emission  increases  with 
density  less  than  with  a  linear  dependence).  As  noted  by  Gadamer  (Ref .30),  at 
approximately  5  mm  Hg.  pressure  for  room  temperature  air,  the  quenching  de¬ 
excitation  is  very  severe.  The  region  of  acceptable  linearity  extends  only 
up  to  a  few  hundred  microns  pressure  at  room  temperature . 

A  detailed  analysis  of  the  beam-excited  emission  in  nitrogen 
has  been  presented  by  Muntz  (Ref.  l4)  and  will  only  be  outlined  here.  The 
excitation  -  emission  path  is  indicated  schematically  as  follows: 
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Internuclear  Distance 

The  detailed  potential  energy  curves  ean  he  obtained  from  .Gilmore  ;(Ref,.  31')  * 
The  energetic  electrons  excite  and  ionize  the  ground  state  nitrogen  molecule 
(H^1-^),  up  to  the  doublet  excited  at  ate  of  the  molecular  ion  -  The 

subsequent  emission  to  the  ground  state  ion  3C^£)  then  comprises  the  first 
negative  system.  It  is  the  D-0  vibration  -  rotation  hand  that  is  used  in  the 
present  investigation,  having  a  handhead  at  39i*tA°- 

Dor  a  gas  in  thermal  equilibrium,  the  population  in  any 
rotational  level  is  proportional  to  .(following  Serzberg,  Def-  32) 

-  1(1+1)  — | 

%  SX  {23  +  1)  e  ,  (1) 

where  3  is  the  molecular  rotational  constant,  h  is  Blanch'  s  constant,  k  is 
Boltzmann's  constant  and  c  is  the  velocity  of  light-  The  combination  Bhc/k 
has  dimensions  of  °K  and  is  called  the  .characteristic  rotational  temperature 
of  the  molecule,  0p-  Dor  the  nitrogen  molecule,  ~  2-89  °i£.  Since  the 
factor  (2  J+l)  increases  linearly  with  J,  the  number  of  molecules  in  the 
different  rotational  levels  goes  through  a  maximum  as  J  is  increased.  This 
*  maximum  occurs  at 


Thus,  as  the  temperature  decreases,  the  population  maximum  moves  toward  the 
lower  rotation  levels. 
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I&e  variation  of  the  intensity  of  the  lines  in  a  vibration- 
rotation  band  as  a  function  of  *3  is  given  essentially  by  the  thermal  distribu¬ 
tion  of  the  rotational  levels.  In  emission,  the  intensity  depends  on  the  mean 
valts  of  (2  JAl)  for  the  unoer  and  lover  states,  as  indicated. 

Q, 

■  X  -JH-J’-I)  ~ 


1  «{*-(< 


+5^~- 


tL)e 


(3) 


•here  C  is  a  constant,  is  the  rotational  state  sun,  and  J'  denotes  the  upper 
state,  J"  denotes  the  Xci*er  state-  Therefore,  the  intensity  distribution 
resent  las  closely  ths  pcp$iIehi-on  distribution  of  the  rotational  levels,  'a  Jiaxi- 
of  intensity  occurring  at  «ith  Jt-ov  decreasing  as  the  temperature 

•decreases,  shown  schematically 


As  pointed  out  by  Huntz  (Bef.  14),  for  the  1st  negative  system 
of  litregejj  the  electron  sain  ccralisg  in  the  2t>+  3~£  state  is  very  weak, 
and  all  transitions  between  electronic  states  are  governed  by  Eund’s  case  (b). 

If  the  doublets  (i-e.  for  the  B^+  3r£  state)  are  unresolved,  the  transitions 
beccoe  equivalent  to  *£  -*  ‘£  transitions  with  rotational  levels  designated  by 
H,  the  relational  caaatca  number  apart  from  spin,  J  =  XZ|r  (see  Bef.  32).  The 
selection  rule  applying  to  Eumd*  s  case  (b)  for  a  *E  -»*£  transition  is  =  +  1. 
This  selection  rule  results  in  the  formation  of  the  two  familiar  P  ana  R  branches 
in  the  rotation  emission  spectrum. 

For  the  3  branch,  the  wave  length  of  emitted  light  decreases 
with  increasing  Ks  (i.e.  the  lines  are  spaced  toward  the  violet  end  of  the 
spectrum)  spacing  between  lines  increases  with  K’ .  The  P  branch  lines  are 
space i  toward  the  red  until  this  braised  doubles  back  on  itself  forming  a  band 
head.  For  the  C—G  vibration-rotation  band  studied  here,  the  band  origin 
(w*  a  o)  is  at  3909  .4°  (Bef.  33)  s  and  the  band  head  due  to  the  F  branch  folding 
back  lies  at  3914  A°.  The  part  of  the  B  branch  used  in  the  investigation  con¬ 
sists  of  about  21  lines  -  the  reason  being  that  the  rotational  band  bead  of 
the  1-1  vibration  transition  occurs  at  snout  3^84  A°,  onus  the  tonal  unobstruct¬ 
ed  wave  length  interval  involved  in  the  B  branch  has  a  span  of  about  25  A°. 

Tills  is  shown  schematically  in  the  sketch  below. 
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If  the  assumption  is  made  that  the  ground  state  nitrogen  mole¬ 
cule  before  excitation  by  beam-electrons  has  a  thermal  equilibrium  distribution 
of  rotational  energy  (see  Eq.  l)  at  a  rotational  temperature  TR,  and  that  the 
electronic  transitions  are  governed  by  optical  selection  rules,  it  is  possible 
to  predict  the  line  intensities  in  the  1st  negative  system.  Muntz's  (Ref .14) 
relation  for  this,  valid  for  temperatures  below  about  800°K,  (to  preclude  any 
appreciable  population  in  all  but  the  lowest  vibrational  level)  is  given  as 


(VAjV*  (K'+K”*1)  X4 


[8] 


V  4  -K'(K'+l)eRl/TR 
■or—  e 


Here,  XR  is  a  constant,  one  prime  indicates  the  upper  electronic  state 
(il>+  ),  two  primes  indicate  the  lower  electronic  state  (lf2+  X2Z  ),  and 

the  subscript  one  indicates  the  ground  state  molecule  (N2  X'E  ).  The  factor 
[G]  ,  which  includes  the  Hbnl-London  rotational  transition  probabilities  is 


-2(K’+1) 


(K'+l)e 

LuJ  - 
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+  K'e 
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2  K'  +  1 


It  is  seen  that  the  factor  [G]  is  itself  a  function  of  Tp. 

From  the  measured  rotational  line  intensities,  the  slope  tech¬ 
nique  (Ref.  32)  can  be  used  to  obtain  the  value  of  Tr.  This  is  seen  by  re¬ 
writing  Eq.  (4)  as 


log 


K'(K’+l)  +  const. 


(6) 


Plotting  the  left  side  versus  K' (K'+l)  will  yield  a  straight  line,  with  slope 
-  8p-j/Tp  ,  from  which  TR  =  Sp^/ slope .  This  technique  is  often  used  to  obtain 
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temperatures  in  high  temperature,  self- radiating  gas  flows  (Ref.  32),  and  in 
the  electron  beam  experiments,  the  energetic  electrons  are  vised  to  stimulate 
emission  in  a  non-radiating  gas.  It  should  be  noted  that  the  procedure  is 
iterative,  in  that  a  guessed  value  of  Tr  has  to  be  used  to  compute  [G]  ,  before 
a  Tr  value  can  be  obtained  from  Eq.  6.  This  value  is  then  vised  to  compute  [G] 
again,  etc.  For  temperatures  not  far  different  from  room  temperature  (i.e. 

Tr  >  150°K)  one  iteration  usually  suffices,  since  [G]  is  not  a  strong-  function 
of  Tr  in  this  range,  and  only  a  small  error  will  arise  if  only  the  first  guess 
on  [G]  is  used.  However,  at  the  low  temperatures  reached  in  free  jets,  [G] 
becomes  a  strong  function  of  Tr  (i.e.  as  Tr  approaches  8r  in  the  exponentials), 
and  a  complete  iteration  is  necessary.  Values  of  log  [G]  (  v/vQy  are  given 
by  Muntz  (Ref.  14)  for  Tjj  >  75  °K,  and  this  compilation  has  been  extended  for 
the  present  study  to  Tr  =  8  °K.  Tiese  values  are  given  in  Table  I  (where  the 
factor  (v/vQ)^  is  ~  1,  due  to  the  small  wave  length  range  involved) . 


The  reduction  of  the  spectral  data  obtained  in  the  experimental 
program  went  as  follows: 

the  heights  (intensity  of  the  rotational  lines) were  measured  with  a 
scale,  then 

[log  (height)-log  (K'+lC'+l)  -  log  (G) ] 
was  plotted  versus  K* (K'+l).  The  slope  of  the  resulting  straight 
line  enabled  the  rotational  temperature  to  be  computed. 

At  low  temperatures,  care  would  be  taken  to  be  sure  that  the  Tr  guesses  (to 
compute  [G]  )  would  be  above  and  below  the  actual  and  a  linear  extrapola¬ 
tion  between  the  two  was  used  to  obtain  the  final  value  of  Tg  (,i.e.  guesses 
of  Tr  =  10  °K  and  TR  =  15  °K,  were  used  to  obtain  the  actual  TR  value  of  13  °K, 
for  example) . 


As  can  be  seen  from  Eqs.  1  and  2,  the  population  of  the  higher 
rotational  levels  falls  off  rapidly  as  the  temperature  is  decreased.  The 
position  of  maximum  intensity,  for  example,  shifts  from  K'  ~  7  at  300  °K,  to 
about  K'  *  2  at  20  °K.  Thus,  the  number  of  intense  lines  obtained  in  a 
spectra  decreases  with  temperature,  from  21  lines  at  300  K,  to  approximately 
7  lines  at  20  °K,  for  example. 

The  technique  described  above  has  been  used  to  obtain  rotational 
temperatures  in  free  jets  and  shock  waves.  The  various  comments  concerning 
position  of  the  maximum  intensity  line  and  number  of  lines  obtained  in  a 
spectrum,  etc.,  can  be  more  fully  appreciated  when  examining  the  spectral  data 
presented  in  Section  IV. 


\ 
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III.  EXPERIMENTAL  APPARATUS 


3.1  Lov;  Density  Wind  Tunnel 

All  experiments  were  performed  in  the  UTIAS  low  density  wind 
tunnel,  a  facility  that  has  been  in  operation  for  a  number  of  years.  The 
original  design  and  operation  study  of  this  facility  was  reported  by  Enkenhus 
(Ref.  34)  in  1957 }  and  since  then,  the  capability  of  the  facility  has  been  in¬ 
creased  with  additional  instrumentation.  An  example  of  this  is  the  electron 
gun  chamber  mounted  on  the  top  of  the  test  section  on  a  traversing  assembly. 
This  allows  a  beam  of  energetic  electrons  to  be  directed  vertically  downward 
in  the  test  section,  with  this  beam  capable  of  traversing  in  a  square  of 
approximately  six  inches  a  side  (i.e.  six  inches  downstream  and  three  inches 
radially  outward  from  the  tunnel  centerline) . 

The  tunnel  itself  is  a  continuous  running  facility  having  a 
primary  pumping  system  of  nine  large  booster  type  oil  diffusion  pumps.  Each 
of  these  boosters  has  a  pumping  capacity  of  approximately  1,000  liters  per 
second  at  tunnel  test  section  pressures  from  1  to  60  microns  Hg.  The  booster 
pumps  are  connected  to  a  second  stage  of  pumping  consisting  of  two  large 
positive  displacement  mechanical  pumps  of  48?  cubic  feet  per  minute  capacity. 
The  test  section  is  a  cylinder  44  inches  in  diameter,  with  its  axis  at  right 
angles  to  the  rest  of  the  tunnel.  A  schematic  of  the  tunnel  can  be  seen  in 
Fig.  4.  A  remotely  controlled  traversing  assembly  inside  of  the  test  section 
was  used  in  the  present  experiments  to  support  and  align  the  shock  holder  in 
the  shock  wave  studies.  The  operating  conditions  for  the  free  jet  experiments 
to  be  discussed  later  covered  the  following  ranges: 

stagnation  pressures  (Po)  from  .5  to  250  mm  Hg 

and  test  section  or  anibient  pressures  from  about  10  to  30  microns  Hg. 


3.2  Electron  Gun  Chamber 


The  traversing  electron  gun  assembly  (Ref.  35)  mounted  atop 
the  tunnel  test  section  is  the  primary  diagnostic  toe]  in  the  present  experi¬ 
ment.  For  a  schematic  representation  of  the  gun  chamber  and  vacuum  system, 
see  F.ig.  5* 


A  Phillips  TV-tube  is  mounted  in  the  inner  of  the  two  concentric 
tubes  shown  in  the  figure.  An  oil  vapor  diffusion  pump  (Edwards  model  203  B, 
rated  at  50  liters  per  second  with  a  water  baffle)  maintains  the  pressure  in 
this  tube  between  2  x  10"°  and  5  x  10"-5  mm  Hg.  An  Edwards  liquid  air  cold 
tr?"  is  placed  between  the  water  baffle  and  the  gun  chamber.  The  outer  tube 
shown  in  Fig.  5  forms  an  intermediate  pressure  chamber  separating  the  inner  gun 
chamber  from  the  test  section  of  the  low  density  tunnel.  A  similar  diffusion 
pump,  water  baffle  and  cold  trap  system  is  used  to  maintain  the  pressure  in 
this  chamber  at  10“4  to  10"3  mm  Hg,  while  the  test  section  pressure  may  be 
typically  in  the  range  from  1  -  100  microns  Hg.  The  inner  and  outer  chambers 
can  be  isolated  from  the  diffusion  pumps  by  gravity-operated  butterfly  valves 
in  the  event  of  a  power  failure. 

The  electron  beam  passes  through  a  2  mm  diameter  nozzle,  which 
is  12.5  mm  long,  into  the  intermediate  chamber,  and  from  there  the  beam  exits 
into  the  test  section  through  a  15  mm  long  nozzle  with  a  1  mm  diameter  hole. 
Careful  alignment  and  focusing  of  the  beam  is  necessary  to  allow  it  to  pass 
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through  both  nozzles.  This  alignment  is  achieved  with -a  set  of  eight  electro¬ 
static  deflection  plates  within  the  inner  chamber .  Independent  controls  are 
provided  for  shifting  or  tilting  the  beam  in  two  perpendicular  planes. 

The  cathode  potential,  is  supplied  by. an  B5E  high-voltage  regulated 
DO  power  supply  and  may  be  varied  between  0  and  30  kilovolts  (for  all  experi¬ 
ments  described  in  this  report,  the  cathode  potential  was  17. 5  kilovolts), 
and  the  beam  current  was  nominally  held  between  100  and  200  microamperes.  The 
anode  is  kept  at  ground  potential  while  the  grid  and  lens  voltages  are  obtained 
from  a  voltage  divider  circuit.  A  separate,  regulated  400  volt  DC  power  supply 
furnishes  the  deflection  plate  voltages.  A  schematic  of  tbs  gun  and  associated 
circuitry  is  shown  in  Fig.  6.  , 

3  <-3  The  Spectrometer 

The  success  of  an  experiment  of  this  type  depends  upon  the 
availability  of  a  sensitive  (required  because  of  the  low  light  output  at  low 
gas  densities),  high  dispersion  spectrograph  (enough  dispersion  to  separate 
all  of  the  lines  in  the  rotational  spectrum,  a  total  range  of  about  25-  %)  - 
The  instrument  was  originally  designed  by  Marsden  to  fulfill  the  above  mentioned 
requirements,  and  is  described  in  detail  in  Ref.  16.  It  is  an  FDert  type  scann¬ 
ing  spectrometer,  using  a  Bausch  and  Lomb  replica  grating  102  run  square  with 
1200  rulings  per  am  blazed  for  3500  Ac  in  second  order.  The  spherical  mirror 
has  a  focal  length  of  750  mm.  This,  together  with  the  size  of  the  grating, 
gives  the  instrument  an  optical  f  number  (f  equals  foeal  length  divided  by 
the  effective  grating  diameter)  of  6.6.  The  instrument  was  use_  second 
order,  with  a  dispersion  of  4.7  X  per  mm.  The  slits  are  curved  to  provide 
sharp  resolution  over  the  full  slit  length,  see  Fastie  (Ref.  36).  A  schematic 
of  the  instrument  is  given  in  Fig.  7- 

The  light  passing  through  the  entrance  slit  is  focussed  by  the 
spherical  mirror  into  a  parallel  beam  of  light  incident  on  the  grating. 

Spectrally  resolved  light  from  the  grating  is  then  focused  by  the  mirror  on  the 
plane  of  the  exit  slit.  Rotation  of  the  grating  (via  the  grating  drive  motor 
and  lever  arm)  causes  the  resolved  spectrum  to  move  past  the  exit  slit  and  the 
intensity  of  the  lines  is  measured  by  a  photomultiplier  placed  behind  the  slit. 
For  the  experiments  discussed  in  this  report,  the  entrance  slit  width  was 
0.004  inches  (100  microns)  and  height  was  0.300  inches,  while  the  exit  slit 
width  was  .006  inches.  As  discussed  by  Marsden  (Ref.  16),  with  tnese  slit  di¬ 
mensions,  all  lines  with  K’  >  3  should  be  resolved  with  no  contribution  (or 
rotational  line  overlap)  from  neighboring  lines.  Lines  K’  =  1,  2,  however,  may 
be  affected  by  the  presence  of  its  neighbor  (since  line  spacing  decreases  with 
decreasing  K').  The  experimental  results  presented  in  Section  IV  seem  to  bear 
out  this  conclusion. 

The  photomultiplier  used  to  measure  the  line  intensities  is  an 
uncooled  EMI  type  9502  S,  having  an  S-ll  spectral  response  (a  peak  at  4000  $) . 

The  dark  current  as  recorded  on  a  Keithley  4l4  micro-microameter  was  found  to 
be  1.5  -»2.0  x  10-1°  amperes  at  1500  volts. 

Careful  alignment  of  the  grating  and  slits  was  accomplished  using 
light  sources  of  mercury,  hydrogen  and  nitrogen.  Proper  slit  alignment  is 
critical  for  experiments  at  low  temperatures  where  only  a  few  spectral  lines 
are  recorded.  These  lines  have  low  K'  values  and  small  spacing,  for  example, 

AA  for  K'  =  1  to  K'  =  2  is  only  O.67  £.  Thus,  the  entrance  and  exit  slits 


10 


y^»v^a^r.^«9r,Is?,iwg-^rg?tWMy  gltH-W^Urrrt^w^yiiRvWHWftt'tWtfPBy^SWWWPWTOffWHfflSWi^S^WPBB^W^WWB^roiBHIBWSWj^H^^P^^C^^^ffWl^^ 


must  be  carefully  aligned  to  minimize  the  amount  of  line  overlap  at  low  values 
of  K* .  This  vas  accomplished  by  removing  the  j&otceultiplier  at  the  exit  slit, 
and  visually  watching  various  lines  from  the  light  sources  pass  by  the  slit  as 
the  grating  was  rotated  by  hand  using  the  micrometer  drive,  Any  non-alignment 
of  entrance  and  exit  slits  could  be  observed  by  a  "rolling"  of  the  emission 
line  as  it  passed  by  the  exit  slit.  The  exit  slit  vas  then  rotated  until  the 
emission  line  jumped  sharply  into  and  cut  of  the  slit.  Proper  alignment  seems 
to  have  been  achieved  in  this  manner,  as  can  be  seen  from  the  spectral  data 
shewn  in  Section  IV. 

3-4  Experimental  Setup 

The  complete  experimental  arrangmeni  is  shown  schematically  in 
Fig.  8,  which  shows  a  top  view  of  the  system.  The  fluorescence  from  the  beam 
is  focused  onto  the  entrance  slit  of  the  spectrometer  with  a  K-24  aerial  camera 
lens.  It  should  be  noted  that  the  entrance  slit  is  parallel  to  the  beam  (i.e. 
beam  vertical,  slit  vertical) .  This  focusing  leas,  which  has  a  focal  length  of 
7  inches,  and  an  optical  f  number  of  2.5,  is  attached  on  a  bar  support  which  is 
part  of  the  election  gun  assembly.  Thus,  the  lens  moves  with  the  beam  as  it  is 
traversed  axially  (or  radially)  in  the  free  jet.  The  lens-to-beam  and  lens-tc- 
slit  distances  were  approximately  9  snd  34  inches  respectively,  giving  a  magni¬ 
fication  of  about  4.  Sins,  the  -004  inch  wide  entrance  slit,  when  imaged  on  the 
1  mm  wide  electron  beam  is  only  .001  inch  in  width  and  about  2  m  in  height. 

This  defines  the  "point"  or  elemental  volume  in  the  beam  fluorescence  where  the 
measurements  were  cade. 


The  spectrometer  and  grating  drive  motor  assembly  rested  on  a 
specially  constructed  traversing  table  which  rode  on  ball  bushings-  This  table 
allowed  the  spectrometer  to  "track"  the  beam,  i.e.  to  "track"  the  beam  image 
formed  by  the  focusing  lens,  to  keep  this  image  on  the  spectrometer  entrance 
slit.  The  lead  screws  for  this  table  were  similar  to  those  for  the  electron 
gun  traversing  mechanism  atop  the  tunnel  test  section,  8  Acme  threads  per  inch. 
Graduated  dials  on  the  lead  screws  were  marked  in  thousandths  of  an  inch  for 
both  traversing  assemblies. 

Several  photographic  views  of  the  experimental  apparatus  can  be 
seen  in  Fig.  9(a)  through  9(d)-  The  first  photo  shows  the  electron  gun  and 
associated  vacuum  system  mounted  on  its  traversing  assembly  atop  the  tunnel  zest- 
section.  This  assembly  moves  on  a  sliding  0-ring  seal  preventing  leaks  from 
occurring  in  the  test  section.  Also  seen  in  the  photo  is  the  mechanical  back¬ 
ing  pump  for  the  diffusion  pumps,  also  mounted  atop  the  low  density  tunnel.  A 
view  of  the  spectrometer  and  traversing  table  in  positron  at  the  test  section 
window  is  given  in  Fig.  9(b)*  The  traversing  table  lead  screws  and  ball  bush¬ 
ing  supports  can  be  seen,  along  with  the  grating  drive  motor  mounted  directly 
on  the  table.  Figure  9(c)  shows  a  general  view  of  the  associated  instrumenta¬ 
tion,  including  the  high  voltage  power  supply,  photo-multiplier  power  supply, 
micro-microammeter,  etc.  The  final  photograph  shows  an  internal  view  of  the 
tunnel  test  section,  as  viewed  from  the  opposite  side  of  the  tunnel  from  the 
spectrometer.  The  bottom  of  the  electron  gun  chamber  with  the  automatic  door 
and  the  beam  collector  cup  are  clearly  seen.  In  position  are  the  focusing 
lens  described  above,  and  the  shock  holder  model  mounted  on  the  test  section 
traversing  mechanism.  The  nozzle  seen  in  the  photograph  is  an  orifice  type  with 
an  exit  diameter  of  5  *• 
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The  system  vas  aligned  optically  in  the  following  manner:  the 
traversing  table  vas  properly  aligned  vith  the  tumal  to  assure  parallel  rove- 
meat  vith  the  electron  beam.  The  spec tresseter  vas  counted  os  the  table  and 
leveled.  A  tungsten  lamp  vas  then  used  to  backlight  the  spectrometer  vith  the 
exit  slit  removed.  An  image  of  the  entrance  slit  vas  forced  (via  the  focusing 
lens)cn  a  ground  glass  screen  inside  the  tunnel  test  section.  A  plucb  line  vas 
suspended  from  the  automatic  dcor  at  the  b ottos  of  the  gun  assembly,  in  line 
vizh  the  electron  bear.  Tee  entrance  slit  image  vas  then  focused  on  the  ground 
glass  at  the  plucb  lice,  2nd  the  spectrograph  height  vas  adjusted  to  insure 
correct  focusing  and  alignment  vlzh  the  centerline  of  the  tunnel.  This  vas 
accomplished  by  moving  the  plucb  line  doss  to  the  exit  orifice  of  the  nozzle 
under  study,  and  eyeball  centering  vas  used.  The  gun  assembly  vas  then  tra¬ 
versed  along  the  eenterline  of  the  tunnel,  tracked  by  the  spectrometer,  to 
insure  proper  focusing  for  the  full  length  of  travel.  This  procedure  vas 
followed  vith  the  installation  of  every  nozzle,  and  vas  felt  to  be  accurate  to. 
less  than  .010  inch. 

The  output  of  the  photomultiplier  vas  amplified  by  a  Xeithley 
hlb  micro-microammeter  and  fed  into  a  Tektronix  ylU  storage  oscilloscope. 

During  the  experimental  runs  the  spectrometer  vas  centered  on  the  bean  by 
observing  the  output  cn  the  micro-mcroarmeter  dial.  Using  one  line  of  the 
spectrum,  the  output  vas  maximized  as  the  entrance  slit  image  vas  traversed 
back  and  forth  across  the  beam  fluorescence  (i.e.  traversing  the  spectrometer 
vith  its  traversing  table) .  Since  the  beam  vidth  vas  about  1  rm  total,  and  the 
slit  image  vidth  only  .001  inch,  this  gave  a  rather  broad  plateau  of  intensity 
of  about  .020  inch.  The  maximum  part  of  this  plateau  vas  considered  as  the 
beam  center  and  the  spectrometer  vas  aligned  vith  this.  After  some  experience 
vith  numerous  runs,  this  maximum  vas  found  to  be  trackable  to  -within  about  .005 
inch. 

The  tins  constant,  or  response  speed,  of  the  Keithley  amplifier 
varied  vith  the  scale  range  being  used.  However,  for  most  of  the  experiments 
the  range  vas  about  ICT^  to  10""^  amperes,  giving  a  time  constant  of  about  ,07 
seconds  to  63$  of  final  current.  Several  grating  scanning  speeds  were  employed 
in  the  early  experiments  to  determine  the  optimum  speed  needed  to  obtain 
accurate  spectra,  and  enough  data  points  for  any  given  free  Jet  or  shock 
wave  setup.  The  final  speed  chosen  gave  a  dwell  time  at  the  peak  of  each  line 
of  about  2  seconds,  much  larger  than  the  time  constant  of  the  amplifier.  These 
results  will  be  discussed  in  Section  IV. 

Before  actual  spectral  data  was  taken  during  a  run,  a  photograph 
of  the  flow  was  taken  using  a  Polaroid  camera  mounted  on  the  opposite  side  of 
the  test  section.  In  this  way,  a  visual  examination  of  the  flow  field  could  be 
made  to  insure  proper  flow  conditions,  etc.  This  beam  visualization  technique 
is  described  by  Rothe  (Ref.  37),  and  consists  in  traversing  the  electron  beam 
parallel  to  itself  down  the  centerline  of  the  Jet.  The  camera  shutter  w^.s  left 
open  and  the  flow  conditions  then  are  “painted"  on  fast  Polaroid  film.  Photo¬ 
graphs  of  underexpanded  Jezs  and  shock  waves  are  shown  in  Section  IV. 
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The  experimental  program  consisted  of  a  systematic  survey  of 
under- expanded  nitrogen  jets  issuing  from  a  room  temperature  source,  using  the 
electron  beam  probe  to  obtain  rotational  temperature  distributions.  A  number 
of  sonic  orifices  and  stagnation  pressures  were  used  in  the  investigation  to 
alloy  a  large  variation  in  flow  properties.  Several  shock  waves  ver^  investi¬ 
gated,  covering  the  Mach  number  range  from  about  M  =  k  to  15-  ?or  these  runs,  a 
shock  holder  model  was  inserted  in  the  free  jet  before  the  first  Hach  disc  (Ref. 
35)-  In  addition  to  temperature  measurements,  density  profiles  were  obtained 
for  both  free  jets  and  shock  waves  using  a  photomultiplier  optical  system  de¬ 
signed  by  Rothe  (Refs.  35  s  38).  Ins  results  of  the  measurements  are  given  in 
the  present  Section  of  this  report,  it  being  divided  into  three  separate  sub¬ 
sections  .  A  discussion  of  the  results  and  conclusions  will  be  presented  in 
Section  V. 

b.l  Calibration  Runs 


After  the  spectrometer  was  optically  aligned  with  the  beam-lens 
system  (as  described  in  Section  III) ,  a  number  of  measurements  were  rads  at 
room  temperature  to  give  an  indication  of  the  accuracy  and  scatter  of  the 
rotational  temperature  data.  Room  temperature  bottled  nitrogen  (Linde)  was 
admitted  directly  into  the  tunnel  test  section  through  a  large  diameter  side 
port.  One  oil  booster  pump  was  used  to  maintain  a  quiescent  nitrogen  ambient 
condition  in  the  chamber  at  a  fixed  pressure,  ?m  .  The  temperature  of  the 
nitrogen  in  the  test  section  was  monitored  with  a  copper- const antan  thermocouple, 
with  the  output  read  directly  on  °F  on  a  Brown  Instruments  potentiometer.  The 
electron  beam  was  then  turned  on,  end  a  rotational  spectra  was  obtained  -  the 
thermocouple  output  being  continuously  monitored  during  the  time  interval  re¬ 
quired  for  a  spectral  scan.  As  described  previously,  the  grating  speed  was 
such  that  the  entire  R-branch  spectrum,  was  obtained  (for  this  room  temperature 
condition)  in  approximately  5  minutes.  The  electron  beam  current  was  also 
monitored  during  this  interval  on  a  microampere  meter  measuring  the  collector 
cup  current.  If  the  current  varied  by  more  than  5p  during  this  time,  the  run 
was  discarded  (this  procedure  was  used  throughout  all  experiments) .  A  sweep 
speed  of  about  10  seconds  per  division  was  used  on  the  oscilloscope,  giving 
sufficient  resolution  such  that  all  the  line  heights  could  be  measured  accurately. 
Thus,  about  3  photographs  of  the  scope  face  were  needed  for  a  given  spectral 
trace.  At  the  end  of  the  each  sweep,  the  grating  drive  motor  was  stopped  mo¬ 
mentarily  while  the  scope  was  re-triggered,  then  the  spectral  trace  was  resumed, 
etc.  This  momentary  interruption  in  the  grating  scan  did  not  affect  the  data, 
as  will  be  seen  from  both  the  actual  spectral  photographs  and  the  line  heights 
used  in  the  straight  line  slope  technique. 

These  calibration  runs  were  performed  at  test  section  pressures 
ranging  from  about  10  to  125  microns  Hg.  The  pressures  were  measured  with  a 
cold  trapped  mercury  McLeod  gauge  mounted  outside  of  the  tunnel.  The  run  con¬ 
ditions  are  summarized  in  the  following  table: 
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Figure  10  shows  a  typical  spectrins  obtained  at  room  temperature  conditions. 
Eighteen  rotational  lines  in  the  R  branch  are  seed,  with  the  spacing  between  lines 
increasing  with  increasing  X* ,  and  K*w;3y  being  about  X*  =  7-  Owing  to  the 
nuclear  spin  properties  of  the  nitrogen  molecule,  the  alternate  line  intensities 
are  predicted  to  be  in  the  ratio  of  2.1  on  the  X*  (X’-j-l)  log  plots  at  the  same 
value  of  the  ordinate  (Ref.  32) .  These  half-intensity  lines  correspond  to  even 
values  of  X’  as  shown  in  the  figure.  For  all  data  reduction,  the  heights  of 
these  lines  were  multiplied  by  two  for  use  in  the  slope  technique,  shown,  for 
example,  in  Fig-  10  for  the  same  spectrum  indicating  a  Tg  =  298°X.  It  is 
seen  that  the  points  are  linear  with  X*  (X’ti)  showing  a  Boltzmann  distribution 
to  exist  in  the  ground  state  molecule.  In  addition  the  even  X’  lines  fall  in 
line  with  the  odd  X*  lines,  indicating  an  alternating  intensify  ratio  close  to 
the  theoretical  figure  of  2.  There  was  50 me  question  about  the  value  of  this 
ratio  previously  (Sef.  l6},  but  careful  slit  alignment  aim  mixed  line  overlap 
that  would  have  affected  this  ratio  experimentally,  and  the  factor  of  2  now 
seems  to  be  correctly  measured. 

Figure  11  shows  the  linear  plot  of  log  intensity  vs.  X' (K’-KL) 
used  to  obtain  the  £3  values  for  the  rcom  temperature  calibration  runs.  It  can 
be  seen  from  the  data  that  a  scatter  of  about  t  2-5 dp  exists  in  the  present  data 
at  300°X,  which  is  in  general  agreement-  with  Muntz’s  (Ref.  l4)  results.  There 
is  some  indication  from  the  log  plots  that  the  first  two  points  cay  be  high. 

This  could  be  due  to  the  overlap  of  neighboring  lines  at  the  low  X‘  values 
(this  was  mentioned  in  Section  III  in  the  discussion  of  spectrometer  slit  di¬ 
mensions}.  There  is,  however,  another  source  that  may  contribute  to  the  apparent¬ 
ly  higher  intensity  of  these  lines.  Since  the  P  branch  doubles  back  upon  itself 
(to  fora  the  bandhead  at  3914  2),  it  cay  be  expected  that  a  number  of  the 
lightly  populated  higher  Xs  lines  of  this  branch  may  be  in  the  approximate  wave 
length  region  of  the  low  X'  lines  of  the  R  branch  under  study.  These  would 
contribute  to  the  measured  intensity  of  the  R  branch  lines.  This  was  discussed 
by  Marsden  (Ref.  16)  for  conditions  of  500°K,  where  the  overlapping  P  branch 
lines  may  contribute  20^  of  the  measured  R  intensities  for  the  first  3  lines. 

A  similar  calculation  was  performed  for  300°K,  and  as  shewn  below,  the 
intensity  overlap  at-  this  temperature  is  extremely  small: 
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The  P  branch  intensity  contribution  will  be  even  less  than  the 
above  values  at  the  low  temperatures  encountered  in  free  jets,  since  the  high 
K‘  levels  vill  have  no  appreciable  population.  Examination  of  the  straight 
line  log  plots  for  a  large  number  of  experimental  runs  at  room  temperature, 
indicated  that  the  measured  intensities  of  the  first  two  lines  seemed  to  be 
about  5/e  high.  Since  this  value  is  larger  than  P  branch  contribution  predicted 
above,  it  can  be  taken  as  a  possible  indication  of  the  amount  of  rotational 
line  overlap  due  to  slit  geometry  as  discussed  in  Section  III. 

Individual  rotational  lines  were  scanned  with  a  rapid  sweep  on 
the  oscilloscope  to  obtain  line  profiles.  Figure  12  shows  such  a  trace  for 
the  K’  =  3  line™  The  sweep  speed  is  one  second  per  division,  indicating  a 
dwell  time  at  the  peak  height  of  about  two  seconds  (see  Section  III).  This 
grating  scan  speed  vas  used  throughout  the  experimental  program. 

4 .2  Free  Jet  Studies  . 

A  large  number  of  nitrogen  free  jets  were  investigated  using 
the  electron  beam  probe,  and  both  temperature  and  density  measurements  were 
obtained.  The  stagnation  pressure  and  orifice  exit  diameter  vas  varied  to 
obtain  flow  conditions  where  the  rotational  degree  of  freedom  was  expected  to 
depart  from  equilibrium  in  the  expansion  on  the  basis  of  South's  calculations 
(Ref.  13).  The  parameter  (PQd),  where  PQ  is  the  stagnation  pressure  in  torr, 
and  d  is  the  orifice  exit  diameter  in  mm,  was  varied  frem 

(P  d)  =  15  torr-mm  to  (P  d)  =  480  torr-mm 

U  G 

These  values  correspond  to  Reynolds  numbers  based  on  the  sonic  exit  conditions 
at  the  orifice,  of 


Red*  =  290  to  9,320 

A  large  contoured  nozzle,  4'?. 5  mm  in  diameter  was  used  for  mea¬ 
surements  in  the  very  low  x/d  range.  Several  orifice  type  nozzles  having  much 
smaller  exit  diameters  were  machined,  and  were  attached  to  the  large  contoured 
nozzle  inside  the  test  section.  This  arrangement  can  be  seen  in  Fig.  9(<i), 
which  shews  the  5  ®m  diameter  orifice  in  place.  These  small  nozzles  were 
machined  with  a  very  thin  wall  at  the  end,  thus  the  nozzles  acted  as  orifices, 
as  follows: 
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The  data  presented  by  Asbkenas  and  Sherman  (Bef.  6)  for  orifices 
having  similar  t/d  ratios  indicate  an  effective  orifice  discharge  c  (efficient 
of  about  1  for  values  of  Be^t  greater  than  100.  Thus,,  the  d  value  used  in  the 
x/d  correlation  is  the  actual  exit  diameter  of  the  orifice,,  given  above. 

The  experiments  were  performed  as  follows  (see  Big. 8):  bottled, 
room  temperature  nitrogen  was  admitted  through  flow  metering  valves  into  the 
stagnation  chamber  cf  the  tunnel  after  the  oil  booster  pumps  had.  evacuated  the 
test  section  to  about  1  micron  Hg.  Sufficient  time  was  allowed  for  the  How 
to  stabilize  and  both  the  stagnation  pressure  (P0)  and  test  section  pressure 
(?  )  were  measured  with  the  mercury  McLeod  gauge.  The  electron  beam  was  turned 
on,  and  a  photograph  of  the  flow  field  was  obtained  with  a  Polaroid  camera  as 
described  by  Bothe  (Bef.  37).  The  beam  was  then  traversed  down  the  centerline 
of  the  jet,  being  tracked  by  the  spectrometer  on  its  traversing  table.  At 
numerous  points  along  the  jet  centerline,  rotational  spectra  were  obtained  as 
described  in  Section  4.1.  Again,  the  beam  current  was  continually  monitored 
to  insure  a  variation  of  less  than  5$  in  the  beam  current  during  the  grating 
scan.  This  proved  to  be  an  upper  limit,  since  for  the  majority  of  the  data 
points,  the  beam  current  fluctuations  were  only  about  1%,  and  in  most  cases, 
not  discernable  at  all  on  the  microampere  meter.  Bor  the  larger  values  of 

(PQd),  about  15  data  points  were  taken  along  the  jet  centerline.  In  addition, 
several  radial  temperature  distributions  were  obtained  in  the  jet  by  traversing 
the  electron  gun  radially  outward  from  the  centerline.  Again,  the  beam  was 
tracked  by  the  spectrometer  which  was  traversed  backwards  away  from  the  tunnel 
to  keep  the  beam  focused  on  the  entrance  slit.  The  stagnation  temperature  was 
monitored  by  a  copper-eonstantan  thermocouple  during  the  run.  Since  it  took 
about  3  minutes  to  obtain  a  complete  spectral  scan  near  the  nozzle  orifice 
(i.e.  Tp  «=  200°K  for  example),  and  about  1-1/2  minutes  to  obtain  a  spectral 
scan  further  down  the  jet  at  low  temperatures  (i.e.  about  20  °K  c*r  so),  a  com¬ 
plete  jet  centerline  survey  could  be  completed  in  approximately  1-1/2  to  2  hours, 
including  flow  visualization. 
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The  important  flow  parameters  of  the  nitrogen  free  jet  runs  are 
tabulated  below: 
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where  x^/d  is  the  Theoretical  distance  to  the  first  Mach  disc,  and  x/ d^ajSi 
is  the  position  where  the  last  spectral  data  was  taken  for  each  run. 


The  jets  observed  covered  a  wide  range  of  pressure  ratios  and  in 
all  cases  the  last  spectral  data  point  was  taken  well  before  the  occurrence  of 
the  first  Mach  disc  as  predicted  theoretically  (Ref.  6).  Actually  the  last 
data  point  for  each  run  was  determined  by  the  local  density  in  the  jet,  the 
optical  system  becoming  light  limited  at  low  densities  (i.e.  the  spectral  trace 
became  quite  noisy,  since  the  photomultiplier  output  was  nearing  the  dark 
current  level).  In  general,  this  light  limiting  point  occurred  at  densities 
equivalent  to  a  pressure  of  about  5  microns  iig  at  room  temperature.  It  may  be 
noted  that  the  density  measurements  to  be  discussed  later  indicated  that  the 
first  Mach  disc  is  quite  diffuse  for  low  values  of(PQd) , 

It  is  of  interest  to  examine  the  photographs  of  the  flow  fie  .Id, 
and  these  are  shown  in  Fig.  13  for  various  values  of  (Pod).  Por  small  (PQd) , 
the  entire  flow  field  Is  quite  diffuse  and  the  free  jets  do  not  seem  at  all 
well  defined.  At  the  larger  (P0d)  values  (and  consequently  higher  pressure 
ratios),  the  jets  become  very  sharply  defined  and  the  barrel  shock  quite  clear. 
The  complete  travel  of  the  beam  is  about  six  inches  giving  an  indication  of  the 
first  Mach  cell  size.  The  slight  uneveness  in  the  quality  of  the  photos  is 
due  to  the  fact  that  the  electron  beam  traversing  screw  was  turned  by  hand. 

Examples  of  spectral  records  from  which  the  values  of  rotational 
temperature  were  obtained  are  presented  in  Fig.  l4  for  a  large  number  of  tem¬ 
peratures  (i.e.  positions  along  the  free  jet  centerline,  for  example)  and  indi¬ 
cate  the  quality  of  the  data.  In  all  records,  the  oscilloscope  sweep  speed  is 
about  10  seconds  per  division.  Readily  apparent  in  the  data  is  the  intensity 
maximum  shift  to  lower  K'  values  as  the  temperature  decreases  (Section  II), 
indicating  the  low  population  of  the  higher  rotational  levels.  The  number  of 
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$2  final  rctaafcmal  temperature  data.  is  clotted  in  Fig.  16  as 
flection  of  distance  ivr  tee  Jet  centerline,  x/d_  For  all  values  of  (?._c) , 
the  fr>  values  fall  ec  ace  ?  =  l.i  i sensrcca c  teanerature  carve  at  g^!1  values 
of  x/d.  2a  data  shears  a  gradnal  departure  free  alas  carve,  and  in  tie  low 
(Fr,d)  case,  tee  rotational  teaceratcre  seers  to  freeze  at  a  relatively  high 
value  in  tie  Jet.  2b  scatter  is  seen  to  be  quite  snail  for  been  crifice  dis- 
reters  useu  in  each  (F©d)  curve.  Hess  results  will  be  distressed  in  Section  ¥, 
bet  it  can  be  said  fm  a  cursory  observation  cf  Fig-  1 £>,  teat  as  (?0d)  is  in¬ 
creased,  tic  mcaticnal  tenmeratare  folloas  tbe  isentropic  7  =  1.4  mrre  to 
1  :wer  temperature: ,  as  would  be  eacected  for  a  ssnc-eqailibritas  phenomena. 


For  all  of  tbe  free  Jst  expe rlmenas,  tee  measured.  centerline  ten- 
peratu'w  departed  freer  tbe  isentrsroic  carve  at  different  locations,  depending 
re.  tee  mine  of  (r^d).  2b  reasrrenerts  obtained  at  tbe  lowest  vales  cf  ?0d  = 

It  t.-rr-iao,  for  example.  .depart  at  about  S3  c5£-  At  tbe  bluest  vales  of  ?0d  = 

^53  t xr:  -cel,  ice  data  seems  to  assart  at  approximately  50  °2,  as  can  be  seen 
free  s  ig.  id-  A  comparison  of  tbe  eager! .rental  data  obtained  sear  tbe  orifice 
(i.s.  at  small  x/d)  before  departure  fbcoo  tbe  isentropic  carve,  is  shewt  la 
Fig.  If  veer  a  data  points  from  all  orifices  and  (?-d)  values  are  plotted  vs. 
x/d.  The  agreement  wish  tie  isentropic  prediction  is  quite  good  and  tbe 
s  latter  obtained  (about  ~  35®)  is  ccngarsble  to  tbe  scatter  ©stained  in  tbe 
room  temperature  calibration  rrs  discussed  in  Secticsn  4.1. 

A  typical  radial  temperature  distribution  is  shown  in  Fig.  13  at 
an  x/i  of  6-3  for  (F0d)  =  23©  torr-sx.  fee  temperature  distribution  in  tbe 
center  of  tbe  Jet  is  seen  to  be  quite  flat,  and  tbe  increase  in  rotational 
temperature  in  tbe  barrel  shock  is  very  apparent. 

As  a  check  on  whether  the  diffuse  Mach  disc  (as  shown  in  flew 
photo.-}  at  low  (?^i)  values  may  be  contributing  to  tbe  departure  of  tbe  rotational 
temperature  from  tbs  7  =  1.4  curve,  a  series  of  runs  were  undertaken  in  which 
the  density  oisiributico  in  the  free  Jets  was  obtained.  The  optical  instrument 
used  is  discussed  in  detail  by  Rotte  (Refs.  35 3  33),  and  will  only  be  described 
briefly  here.  Jt  consists  of  a  phot crralti oiler  housing  counted  to  the  electron 
gun  a ssenbuy  (i.e.  mounted  to  the  sane  oar  used  to  hold  the  focusing  lens  in  the 
speetrome-er  experiments),  so  that  it  traverses  with  the  bean.  Fee  entrance 
slit  to  the  housing  was  focused  perpendicular  to  the  bean  (i.e.  bean  vertical, 
slit  norixontal)  so  that  all  of  the  fluorescence  surrounding  the  bean  would  be 
gathered.  Light  fibres  (iron  the  American  Optical  Company)  were  used  to  guide 
the  radiation  fren  the  entrance  slit  to  an  interference  filter  directly  ahead 
of  the  photomultiplier  face.  An  uncooled  E2£L  9502  S  photomultiplier  was  used 
(similar  to  that  in  the  spectrometer) ,  in  conjunction  with  an  interference 
filter  of  109  2  half -width  centered  at  3900  2.  The  instrument  was  calibrated 
in  the  tunnel  test  section  at  room  tempei ature  in  a  quiescent  nitrogen  atnos- 
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there  -  to  obtain  a  curve  of  light  output  versus  nitrogen  density  (pressure) . 

The  photcmul tiplier  current  was  read  on  the  Keithley  micro-si  croanneter.  This 
calibration  data  is  shewn  in  Fig.  19  and  is  seen  to  be  linear  to  at  least  120 
microns  pressure. 

Using  this  instrument,  a  number  of  free  jets  were  surveyed  to  obtain 
venterline  density  distributions  for  the  5  szs  diameter  orifice.  This  data, 
which  is  shewn  in  Fig.  20,  was  obtained  using  the  calibration  curve  of  Fig.  19. 
The  scatter  is  seen  to  be  quite  small,  and  the  density  distribution  generally 
follows  the  iseatrcpic  curve.  For  low  (PQd)  values,  the  density  leaves  this 
isertropic  curve  and  begins  to  increase  well  ahead  of  the  theoretical  Kach  disc 
location.  This  is  an  indication  that  this  Kach  disc  is  not  at  all  well  defined 
(as  sees  in  Fig.  13) -st  lew  densities,  but  becomes  sharper  as  the  density  is 
increased  at  higher  (?Dd)  values. 

By  comparing  Figs-  lo  and  20,  it  is  seen  that  the  last  spectral 
da ,ta  point  was  tales  before  the  density  begins  to  depart  from  the  isentropic 
carve-  In  fact,  the  rotational  temperature s  shew  marled  departures  from 
equilibrium  while  the  density  still  follows  the  isentropic  curve:  (i.e.  the 
censity  drops  off  approximately  as  (r/d)“-  (Hef .  6)  ).  This  camaiisca  was 
carried  out  fer  the  cost  conservative  case,  since  the  2  mm  diameter  orifice, 
having  such  larger  pressure  raiios9  has  Mach  disc  locatic-ns  farther  downstream 
(ever  twice  as  far  downstream  as  the  5  an  nozzle  in  most  cases:  see  table 
for  free  jet  runs).  Thus,  the  Mach  disc  does  not  appear  to  contribute  to  the 
overall  departure  of  the  rotational  temperature  from  the  isentropic  curve. 

A  more  comprehensive  discussion  of  these  results  is  presented  in 

Section  F. 

h-3  Shoe!  wave  Studies 

•taring  the  course  of  the  experimental  program,  a  shod  holder 
was  inserted  in  the  free  jet  ahead  of  the  first  Mach  disc  so  that  shed  struc¬ 
ture  data  could  be  obtained  for  a  wide  range  of  Kach  numbers .  This  shod  holder 
consists  of  a  hollow,  diverging  cons  vita  a  2-inch  diameter  opening  at  the 
front,  see  Sef-  35-  A  sliding  conical  plug  is  used  to  partially  close  the  rear 
opening,  producing  a  standing  shod  wave  ahead  of  the  holder.  The  leading  edge 
of  the  holder  has  two  sets  of  slits,  i/8  inch  wide  by  1/2  inch  long.  One  set 
of  slits  allows  the  electron  beam  to  pass  vertically  down  through  the  shed 
holder,  and  the  other  set  enables  the  beam  fluorescence  to  be  gathered  by  the 
focusing  lens.  The  shod  holder,  can  be  seen  mounted  in  the  test  section  in 
Fig.  9(d) .  Both  rotational  spectra  and  density  measurements  were  obtained 
through  the  shod  waves  at  closely  spaced  intervals  of  x/d.  The  shod  wave 
runs  are  listed  belcrw: 
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As  in  the  free  jet  experiments,  flow  visualization  photos  were  taken  at  the 
start  of  each  rim.  These  ere  shown  in  Fig.  21 .  It  can  he  seen  from  the  photos 
that  the  shock  wave  is  not  plane.  This  is  due  t-o  the  fact  that  the  free  jet 
flow  is  similar  to  flow  from  a  source  near  the  orifice,  (Ref.  6)  and  the  curved 
shock  then  is  normal  to  the  local  flow  streamlines,  being  perpendicular  to  the 
shock  holder  only  near  the  jet  centerline,  where  the  shock  profile  data  were 
obtained. 


Using  the  photomultiplier-interference  filter  assembly  described  in 
the  previous  section,  density  profiles  were  obtained  along  the  jet  centerline 
and  through  the  shock  wave.  These  are  given  in  Fig.  22  for  the  M  =  4  and  11 
shocks.  The  density  drops  off  similar  to  the  isentropic  curve,  tnen  increases 
through  the  shock  front  towards  the  post  shock  value.  Similar  shock  density 
profiles  have  been  reported  by  Sobben  and  Talbot  (Ref.  39)  •  When  the  electron 
beam  was  near  the  shock  holder,  and  at  times  behind  the  leading  edge,  the 
light  gathered  by  the  focusing  lens  v&s  diminished  due  to  the  presence  of  the 
shock  holder.  To  account  for  this,  an  occulation  calibration  was  carried  out 
in  quiescent  nitrogen  (no  shock  wave),  the  fluorescence  being  measured  as  the 
beam  moved  close  to  the  shock  holder.  This  calibration  curve  is  shown  in  Fig. 
23.  The  data  agrees  extremely  well  with  a  point  source  theoretical  curve 
reported  by  Rothe  (Ref .  38) .  This  calibration  factor  was  used  to  obtain  the 
final  density  data  shown  in  Fig.  22. 


With  the  spectrometer  setup  as  previously  described,  rotational 
spectra  were  obtained  along  the  jet  centerline  through  the  shock  front.  The 
log  slope  plots  are  non-linear,  indicating  a  non-Boltzmann  distribution  in  popu¬ 
lation  of  the  rotational  levels  existing  in  the  nitrogen  molecule  before  excita¬ 
tion  by  beam  electrons.  This  effect  seems  to  be  a  direct  function  of  shock 
strength.  The  deviation  from  linearity  is  small  for  the  M  =  4  shock,  but 
exceedingly  strong  for  the  M  -  15  shock.  Robben  and  Talbot  (Ref.  40)  report, 
for  example,  essentially  no  deviation  at  M  =  1.7,  but  a  much  stronger  deviation 
at  M  =  13.  Figure  24  shows  the  log  slope  distribution  through  a  shock  wave 
for  the  three  cases  investigated,  M  =  4,  11  and  15 -  In  all  cases,  the  log  slope 
is  non-linear  in  the  center  of  the  shock  front,  and  becomes  linear  again  in  the 
high  temperature  region  at  the  tail  of  the  shock  wave.  Thus,  a  Boltzmann  dis¬ 
tribution  is  regained  in  the  post  shock  conditions.  This  non-Boltzmann  effect 
would  go  undetected  in  density  measurements  alone. 

Since  the  factor  (G)  is  a  function  of  rotational  temperature,  T5, 
a  difficulty  arises  when  a  non-linear  log  slope  plot  is  obtained  as  to  what 
value  of  i'g  should  be  used  for  the.  calculation  of  (G) .  Since  the  plots  given 
in  Fig.  24  are  meant  to  be  representative  of  the  non-Bo3.tzmann  effect  with  in¬ 
creasing  shock  strength,  only  one  value  of  (G)  (i.e.  one  value  of  assumed  tr) 
was  used  for  each  slope.  To  obtain  actual  temperatures  from  different  slopes 
of  the  same  log  plot,  detailed  iterations  using  various  assumed  values  of  Tg 
are  needed.  This  procedure  will  be  discussed  in  Section  V.  Thus,  for  the  log 
plots  in  Fig.  24,  the  following  values  of  assumed  Tr  (i.e.  to  compute  (G)  ) 
were  used:  (the  circled  points  denote  positions  in  front  of  the  shock,  @  ,  to 
positions  behind  the  shock,  i.e.  ©  ,  see  Fig.  27). 
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As  pointed  out  by  Marsden  (Ref.  16),  there  is  an  approximate  one- 
to-one  relationship  between  the  line  intensities  in  the  rotational  spectrum  and 
the  rotational  energy  distribution.  Thus,  from  the  plots  in  Fig.  24,  it  appears 
that  the  lower  rotational  levels  have  a  population  reflecting  the  temperature 
anead  of  the  shock,  while  the  higher  rotational  levels  are  populated  according 
to  the  temperature  behind  the  shock  wave.  Rcbben  and  Talbot  (Ref.  40) ,  for 
example,  found  it  was  possible  to  represent  their  results  approximately  by  the 
merging  of  two  rotational  distribution  functions  corresponding  to  temperatures 
upstream  and  downstream  of  the  shock  wave. 


V.  DISCUSSION  AND  CONCLUSIONS 


The  rotational  temperature  measurements  obtained* in  underexpanded 
jets  indicate  that  the  electron  beam  fluorescence  probe  can  be{  used  with  some 
degree  of  confidence  as  a  diagnostic  tool  at  low  temperatures.  As  noted  in 
Fig.  17,  for  example,  the  rotational  temperature  values  obtained  experimentally 
agree  with  the  predicted  centerline  isentropic  distribution  down  to  about  50  °K 
for  high  values  of  (Pod) .  For  all  of  the  experiments,  the  measured  rotational 
temperatures  departed  from  equilibrium,  the  degree  of  departure  depending  on 
the  value  of  (Pcd).  However,  as  noted  earlier,  the  trend  of  the  rotational 
temperature  measurements  with  increasing  (PQd)  is  such  that  the  experimental 
results  follow  more  closely  the  isentropic  curve  to  lower  temperatures.  This 
is  to  be  expected,  since  at  higher  densities  (higher  PQd  values  here),  collisions 
become  more  frequent  and  the  rotational  degree  of  freedom  is  more  apt  to  follow 
the  translational  degree  to  lower  temperatures. 

The  appearance  of  the  rotational  temperature  versus  x/d  curves. 

Fig.  16,  (i.e.  slow  departure  of  rotation  from  equilibrium  with  subsequent  tend¬ 
ency  to  freeze)  is  typical  of  vibrational  non-equilibrium  in  expansion  flows, 
for  example,  and  as  noted,  the  trend  with  pressure,  (PQd) ,  is  correct.  In 
addition,  for  each  value  of  (?Qd),  two  different  orifice  diameters  were  used, 
so  that  for  a  given  Reynolds  number,  i.e.  given  (PQd),  different  flow  densities 
would  be  obtained  (a  factor  of  2-1/2  in  the  densities).  As  can  be  seen  from 
the  small  scatter  for  each  (?0d)  value,  the  data  scales  in  a  binary  fashion. 

The  possibility  of  Mach  disc  interference  with  the  free  jet  rotational 
temperature  measurements  should  be  considered  if  measurements  are  made  far 
downstream.  However,  the  density  distributions  shown  in  Fig.  20  seem  to  rule 
out  this  possibility  for  the  present  results.  Even  at  distances  (x/d)  where 
there  are  large  departures  of  rotational  temperature  from  equilibrium,  the 
density  continues  to  follow  the  isentropic  curve,  decreasing  approximately  as 

(x/d)  • 


It  is  of  interest  to  compare  the  results  of  the  present  investi¬ 
gation  with  the  only  other  published  rotational  temperature  measurements  in 
under- expanded  jets,  those  of  Robben  and  Talbot  (Ref.  21).  Their  experimental 
setup  was  somewhat  similar  to  that  shown  in  Fig.  8,  with  their  electron  beam 
vertical,  but  their  spectrometer  slit  horizontal  (i.e.  across  the  beam).  A 
number  of  measurements  were  made  at  several  values  of  (?Qd),  but  none  were 
carried  out  as  extensively  for  a  given  (PQd)  value  as  those  reported  here.  In 
addition,  their  electron  beam  probe  was  used  at  a  high  beam  current  (approxi¬ 
mately  5  millianperes)  as  compared  to  about  100  microamps  in  the  present  in¬ 
vestigation.  Figure  25  shows  the  University  of  California  Tp  results  compared 
with  those  obtained  at  Toronto.  The  comparison  for  similar  (P0d)  values  is 
quite  good  for  the  two  different  investigations,  but  in  general  it  can  be 
said  their  measurements  somewhat  are  higher  than  those  reported  here.  Robben 
and  Talbot  report  high  Tp  data  for  all  of  their  experiments,  even  those  per¬ 
formed  at  subsonic  Mach  numbers.  The  reason  for  this  is  not  clearly  brought  out 
in  the  reference.  It  is  seen  that  all  of  their  data  lies  fairly  close  to  the 
isentropic  curve,  and  they  have  concluded  that  the  electron  beam  technique  is 
in  error  at  low  temperatures  since  their  measured  rotational  temperatures  are 
above  the  predicted  isentropic  values.  For  example,  when  the  isentropic  tempera¬ 
ture  is  predicted  as  10  °K,  their  Tp  value  obtained  from  spectral  data  is 
about  13  °K,  or  30$  high.  This  is  then  referred  to  as  a  calibration  point, 
i.e.  TR  data  must  be  multiplied  by  0.70  to  obtain  correct  values  in  this 
temperature  range.  However,  the  present  data  shows  that  the  measured  rotational 
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temperature  can  have  numerous  values  depending  on  the  value  of  (?Qd) ,  when 
the  isentropic  prediction  is  10  °K. 

Robben  and  Talbot  (Ref.  21)  do  point  out  one  possible  source  of 
error  associated  vitn  the  electron  bean  technique,  this  is  the  effect  of 
secondary  electrons  on  the  measured  rotational  temperature  (i.e.  the  electrons 
stripped  free  lip  to  yield  2!p~  by  collisions  with  the  primary  bean  electrons). 
These  secondary  electrons  will ,  in  general,  have  much  lower  energies  than 
those  in  the  primary  bean  (3ef.l4),  but  their  actual  energy  distribution  and 
direction  are  unknown. 

Excitation  to  the  excited  ion  state  by  secondary  electrons  may 
arise  free  direct  excitation  of  the  ground  state  molecule  (similar  to  primary 
excitation)  or  excitation  of  the  ground  state  ion.  The  potential  energy 
curves  schematically  diagrammed  in  Section  II  (and  given  in  detail  in  Ref.  3i) 
show  that  electron  energies  in  excess  of  18.7  eV  are  needed  to  excite  the 
ground  state  nitrogen  molecule  (&£  Xl£)  to  the  excited  ion  state  (Hp+  B^E) 
while  electrons  with  energies  in  excess  of  3-1  eV  are  needed  to  excite  the 
ground  state  ion  (Rg”  X^£)  to  the  excited  ion  state.  Thus,  as  pointed  out  by 
Muntz  (Ref.  14)  for  single  collisions,  secondaries  with  energies  below  3-1  eV 
cannot  produce  any  excitation.  In  addition,  Muntz  calculated  that  the  second¬ 
ary  excitation  of  ground  state  ions  (i.e.  for  electrons  with  energies  in 
excess  of  3-1  eV)  may  be  neglected.  For  electrons  with  energies  in  excess  cf 
l8-7  eV,  say  in  the  100  eV  range,  Muntz  indicated  than  there  could  be  a  large 
contribution  to  the  observed  emission  owing  to  direct  excitation  of  the 
ground  state  molecule  by  these  electrons.  If  this  excitation  is  different  in 
nature  than  the  primary  excitation  (i.e.  as  suggested  by  Robben  and  Talbot 
(Ref.  21)  that  these  excitations  say  not  obey  the  optical  selection  rules), 
an  error  could  be  associated  with  the  rotational  temperature  measurements. 

The  fluorescence  excited  by  secondary  electrode  is  confined  to  a 
halo  surrounding  the  beam  of  primary  electrons.  Experiments  .erformed  by 
Muntz,  Abel  and  Maguire  (Ref.  17)  in  helium  indicate  the  presence  of  a  large 
number  of  secondary  electrons,  which  appear  to  be  concentrated  near  the  center 
of  the  beam.  As  pointed  out  in  Ref.  17,  the  contribution  of  secondary  electrons 
to  the  emission  excited  by  an  electron  beam  must  be  determined  as  a  function  of 
radial  distance  from  the  beam  axis  for  each  gas.  Thus,  the  magnitude  of  a 
possible  error  in  the  measurements  may  depend  on  how  much  of  the  halo  is  ob-  _ 
served.  The  University  of  California  experiments,  with  the  slit  perpendicualr 
to  the  beam,  would  gather  more  of  this  halo  emission.  Robben  and  Talbot 
report  that  their  measurements  seemed  to  show  smaller  errors  when  the  spectro¬ 
meter  slit  length  was  reduced  (i.e.  less  halo  observed).  For  the  experimental 
data  reported  here,  this  "effective  heating"  due  to  possible  secondary  electrons 
is  not  apparent  since  the  experimental  temperatures  agree  with  the  isentropic 
curve  down  to  low  temperatures .  This  could  be  due  to  the  fact  that  the  spectro¬ 
meter  slit  is  parallel  to  the  beam  (i.e.  views  a  narrow  slice  along  the  beam 
axis)  and  only  a  small  portion  of  the  halo  within  the  optical  cone  is  observed. 

A  study  of  the  excitation  processes  due  to  secondary  electrons  is  needed  to 
more  fully  understand  the  electron  beam  technique  at  low  temperatures  and 
densities. 


There  is,  however,  another  trend  that  can  be  observed  in  the  pre¬ 
sent  rotational  temperature  data  obtained  in  free  jets.  As  seen  in  Fig.  15, 
the  log  slope  plots  of  rotational  line  intensities  are  linear,  using  the 
majority  of  observed  rotational  lines  at  each  data  point.  However,  the  higher 
K'  rotational  lines  (these  lightly  populated)  always  indicated  a  higher  intensity 
than  was  consistent  with  this  straight  line.  This  effect  was  also  noted  by 
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Robben  and  Talbot  (Ref.  21)  in  a  number  of  their  published  log  plots.  For  a 
given  rotational  temperature  as  obtained  from  the  straight  line  slope,  it  is 
always  the  higher  rotational  lines  that  show  this  apparent  discrepancy  of  over 
population.  Careful  examination  of  all  of  the  free  jet  spectral  data  obtained 
in  this  investigation  indicated  that  the  higher  rotational  levels  seemed  to 
depart  from  a  Boltzmann  equilibrium  with  the  lower  rotational  levels  in  the 
exapnsion.  This  departure^can  be  depicted  in  the  following  way:  define  a 
rotational  "temperature"  Tr^i  based  on  the  slope  for  each  rotational  line  X* 
referred  to  the  K*  =  1  line  as  shown  schematically 


For  a  linearjplot,  (i.e.  Boltzmann  distribution  in  population  of  rotational 
levels),  all  Tr^,  are  equal  and  equal  to  the  rotational  temperature  Tr.  Thus, 
Trk,  is  an  indication  of  the  population  of  a  particular  level  referred  to  the 
lowest  rotational  level.  These  "tempei atures"  have  been  obtained  from  the 
spectral  log  plots  (similar  to  those  shown  in  Fig.  15)  and  plotted  in  Fig.  2 6. 
Again,  where  possible  for  each  (P0d)  value,  results  from  two  different  orifices 
have  been  plotted,  showing  the  same  binary  scaling  as  before.  A  general  trend 

for  each  (P0d)  curve  seems  to  be  a  gradual  peeling  away  of  the  Tr^,  values  for 

the  higher  levels  during  the  expansion^  The  data  indicates  for  example ,  that 
Tr5  peels  away  at  about  17°K,  and  the  Tr7  data  peels  away  at  approximately 
4o°K.  Above  this  peeling  away  point  for 'each  K' ,  these  rotational  lines  lie 
on  the  linear  portion  of  the  log  plot.  It  can  be  seen  for  example,  that  for 
Pcd  =  15  torr-mm,  the  K'  =  7,%, 3  points  never  depart,  since  the  data  does  not 
go  below  1jO°K  (seemingly  the  Tr7  departure  point).  Again  for  Pcd  =  80  torr~mm 

the  K'  =  5,3  curves  never  depart  since  the  data  is  higher  than  17°K  (seemingly 

the  Tr^  departure  point) . 

It  may  be  surmised  that  there  is  no  reason  to  expect  that  there 
would  not  be  a  departure  from  a  Boltzmann  distribution  in  a  rapid  expansion. 

A  greater  amount  of  theoretical  and  experimental  work  has  been  devoted  to  this 
problem  for  the  case  of  vibrationally  relaxing  flows.  The  calculations  of 
Treanor  (Ref.  10),  for  example,  indicate  the  possibility  of  a  non-Bolt. zmann 
population  distribution  in  vibrational  levels  during  expansion.  However, 
experimental  results  of  vibrational  relaxation,  such  as  those  of  Hurle,  Russo 
and  Hall  (Ref.  ll),  have  as  yet  not  confirmed  this.  In  addition,  the  electron 
beam  measurements  of  vibrational  relaxation  in  arc-tunnel  nozzles  reported  by 
Sebacher  and  Duckett  (Ref.  19)  and  Petrie  (Ref.  20)  seem  to  agree  with  the 
results  obtained  using  a  line-reversal  technique  employed  in  Ref.  11.  Further- 


more,  in  the  electron  team  measurements,  a  Boltzmann  distribution  in  vibration 
is  assumed  in  order  to  relate  the  relative  intensities  of  vibrational  band 
systems  to  the  vibrational  temperature - 

This  aspect  of  the  rotational  temperature  measurements  leads  to 
the  possibility  that  the  apparent  rotational  non-Boltzmann  distribution  indicated 
during  the  expansion  process  may  he  produced  by  the  electron  beam  itself,  per¬ 
haps  by  secondary  electrons  in  some  way.  The  author  has  found  no  studies  or 
computations  concerning  the  relaxation  of  individual  levels.  This  would, 

of  course,  require  the  individual  transition  probabilities  of  each  level  to 
be  known,  so  that  a  separate  rate  equation  could  be  applied  to  each  level  during 
a  relaxation  process. 

There  is  a  possible  way  to  determine  experimentally  whether  or 
not  this  departure  from  a  Boltzmann  distribution  is  caused  by  the  beam  or  is 
due  to  the  expansion  process.  This  would  constitute  shooting  an  electron  beam 
into  a  chamber  of  stagnant  low  temperature  nitrogen,  in  equilibrium,,  with  its 
condensate  at  a  very  low  temperature  ( say  30°K)  and  obtaining  "a  rotational 
spectra.  The  log  slope  plot  of  the  line  intensities  should  reveal  any  non- 
Boltzmann  effects  due  to  secondary  electrons.  In  addition,  the  rotational 
temperature  obtained  can  be  compared  to  the  known  temperature  as  a  calibration 
check  on  the  accuracy  of  this  technique  at  low  temperatures.  An  experiment  of 
this  type  has  been  reported  by  Williams  (Ref.  39)  where  an  electron  beam  was 
used  in  conjunction  with  a  liquid  nitrogen  cooled  chamber.  Temperatures  as 
low  as  80°K  were  obtained,  with  the  beam  fluorescence  being  monitored  by  a  set 
of  narrow  pass-band  filters  within  a  single-vibration-rotation  bana,  as 
originally  described  by  Muntz  and  Abe}-  (Ref  ^  4o) .  Unfortunately,  the  reduction 
of  data  with  this  technique  utilizes  an  assumed  Boltzmann  distribution  in 
rotation.  However,  it  has  been  reported  to  the  author  that  additional  experi¬ 
ments  using  a  high  dispersion  spectrograph  to  monitor  the  beam  fluorescence 
are  being  planned. 

The  shock  wave  profiles  presented  in  Fig.  24  can  also  be  repre¬ 
sented  by  a  Trk.  for  each  rotational  line.  In  order  to  make  these  more  de¬ 
tailed  calculations  of  "temperatures"  for  a  non-linear  log  plot,  iterations 
with  (G),  using  various  assumed  values  of  Tr  were  employed.  For  the  Tr^-i 
values  representing  the  higher  rotational  levels  (i.e.  K’>  9)>  the  Tg  values 
given  in  Section  IV  were  used  to  compute  (G).  For  the  lower  levels,  complete 
^iterations  using  appropriate  Tr  guesses  yielded  more  accurate  results  for  these 
TRgt  values.  The  results  are  shown  in  Fig.  27.  The  spread  of  "temperatures" 

Trk>  is  small  for  M  =  4,  since  the  departure  from  linearity  in  the  log  plot  is 
small.  However,  as  the  Mach  number  increases,  the  spreading  is  readily  apparent 
due  to  the  highly  non-linear  log  plots.  For  all  three  cases,  a  Boltzmann  dis¬ 
tribution  (i.e.  all  Trk*  are  equal)  is  regained  in  tail  of  the  shock  wave. 

Since  the  rotational  energy  distribution  in  the  shock  waves  are 
markedly  non-Boltzmann,  one  may  suspect  that  this  effect  could  be  felt  upstream 
of  the  shock  and  perturb  the  rotational  distribution  in  the  jet,  possibly  lead¬ 
ing  to  the  apparent  non-Boltzmann  effect  discussed  earlier.  Robben  and  Talbot 
(Ref.  4l),  for  example,  have  shown  from  their  results  that  the  average  local 
temperature  was  found  to  precede  the  density  profile  through  the  shock  wave  by 
an  amount  which  increased  with  Mach  number.  However,  the  presence  of  a  shock 
wave  in  the  flow  did  not  affect  this  apparent  non-Boltzmann  distribution  in  the 
jet,  as  can  be  seen  from  Fig.  27(c).  For  the  given  (P0d)  flow  condition,  the 
position  of  the  Mach  disc  is  predicted  to  be  at  xj^/d  =  28  and  xjj/d  =  for 
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P  =  torr  and  PQ  =  240  torr  respectively.  She  tzo  sets  of  results,  however, 
agree  very  veil.  In  addition,  the  presence  pi  the  H  =  Ip  shock  caused  by  the 
shock  holder  at  x/d  =  35  ior  the  ?c  =  240  torr  case  did  not  change  the  Tg;r« 
distribution.  Within  the  experimental  scatter,  the  Tgg.  curves  with  the  shock 
vave  present  agree  with  the  pure  jet  results.  Shis  seems  to  indicate  that  the 
presence  of  a  shock  (either  the  Mach  disc  or  a  shock  caused  by  the  shock  holder) 
did  not  influence  this  annarent  non-Boltzmann  distribution  in  the  expansion. 


Another  interesting  aspect  of  the  shock  vave  structure  is  the 
density  profile  data  presented  in  Fig.  22,  which  indicates  that  the  bean 
fluorescence  probe  can  be  used  to  measure  shock  vave  thickness.  So  compare 
the  shock  thickness  obtained  for  the  M  =  3-9  and  M  =  11.5  shock  vaves  presented 
in  Fig  22,  reference  will  be  made  to  the  procedure  of  Hobbes,  and  Talbot  (Bef- 
42) .  A  maximum  slope  thickness  was  obtained  by  drawing  a  straight  line 
through  the  maximum  slope  of  the  density  profile,  and.  measuring  the  distance 
Between  the  intersections  of  this  line  with  the  extrapolations  of  the  tip stream 
and.  downstream  flow  conditions,  i.e.: 


Shis  length  was  then  normalized  with  the  length  =  p(T*)/pu,  where  T*  is 
the  sonic  temperature.  As  pointed  out  by  RoYoen  end  Talbot,  using  S*  allows  a 
more  accurate  computation  of  viscosity  than  when  the  local  mean  free  path  is 
used,  which  requires  a  viscosity  calculation  at  extremely  low  temperatures, 
wher  ■  it  is  cot  well  known.  For  the  computation  of  L*,  the  isentropic  Mach 
number  at  the  shock  holder  leading  edge  was  used  as  discussed  in  Section  3*3- 

Using  p(S*)  -  p  (250  °K)  =  15-45  x  10"^  poise;  li*  was  computed  as: 


L*  cm 


d  nun 
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Referring  to  Fig.  22,  the  maximum  slope  thickness  was  obtained  as: 


At  'Che  lower  Mach  number,  the  thickness  seems  to  agree  veil  with  the  othisr 
data.  However,  at  M  =  11.5,  the  present  L*/^  data  is  about  l/2  of  that  obtain 
ed  fay  Hofafaen  and  Talbot,  but  very  close  to  that  of  Cemac  (Ref.  27)  normalized 
with  L* .  As  pointed  out  in  Ref.  42,  the  large  discrepancy  between  Camac’s 
data  and  that  of  Rofaben  and  Talbot  is  not  clear.  The  indication  here,  however, 
is  that  the  present  data  follows  that  obtained  by  Caxaac,  who  measured  the 
scattered  electrons  from  an  electron  beam  probe  in  a  shock  tube. 


One  other  possible  source  of  error  in  the  free  jet  temperature 
measurements  could  be  condensation  effects  at  low  temperatures  -  that  is, 
molecules  coming  together  to  form  agglomerates  consisting  of  large  numbers  of 
molecules .  These  may  be  expected  to  form,  since  the  conditions  of  low  tem¬ 
perature  and  pressure  that  exist  in  a  free  jet  expanding  from  room  temperature 
are  well  below  the  nitrogen  triple  point  (i.e.  68  °K  and  96.4  am  Hg).  However, 
the  pressures  are  so  low  in  a  free  jet,  that  the  expansion  occurs  such  that 
nitrogen  remains  a  gas  to  a  much  lower  temperature  than  the  triple  point 
temperature  indicated  in  the  sketch^ 


The  vapor  pressure  of  nitrogen  is  a  strong  function  of  temperature 
at  these  lav  temperatures, as  shovn  in  Fig*  28  (taken  from  Ref.  43)  -  Using 
this  figure,  and  the  isentropic  {X  -  1.4)  temperature  and  density  relations 
for  a  free  jet,  possible  condensation  points  were  obtained  for  various  values 
of  stagnation  pressure.  These  are  shovn  in  Fig.  29- 

The  condensation  process  involves  three  body  collisions  in  the 
formative  stages  (i.e.  the  third  particle  carries  away  the  heat  of  sublimation, 
thus  raising  the  gas  temperature  slightly),  so  will  be  a  strong  function  of 
density  (i.e.  a  n^) .  Thus,  higher  stagnation  pressures  will  foster  the  onset 
of  condensation-  Greene  and  Milne  (Ref.  44),  for  example,  report  observations 
of  polymeric  species  for  a  large  number  of  gases,  including  nitrogen,  in  a 
supersonic  molecular  beam  system  expanding  from  room  temperature.  Small  effects 
were  noted  at  1  atm.  stagnation  pressure,,  with  more  pronounced  polymeric  action 
at  5  atm.  pressure.  Additional  skimmer  experiments,  performed  by  Bier  and 
Hagena  (Ref.  45)  indicate  the  onset  of  condensation  at  higher  stagnation  pressures. 
TWitb  a  nozzle  diameter  of  .05  ms,  an  upper  limiting  value  of  (PQd)  =  500  torr-mm 
for  nitrogen  was  found  to  yield  non-condensed  beams.  This,  however,  means  a 
stagnation  pressure  (P0)  of  about  13  atm.  Asbkenas  and  Sherman  (Ref.  6)  re¬ 
port  that  for  stagnation  pressures  ip  to  1  atm.  no  condensation  effect  has 
been  seen  on  inpact  probe  measurements-  A  compilation  of  condensation  mea¬ 
surements  obtained  in  wind  tunnels  and  arc  heated  nozzles  is  given  by  Cassanova 
and  Stephenson  (Sef.  46),  which  indicates  a  high  degree  of  super-saturation  at 
low  temperatures . 

Since  the  condensation  process  is  a  strong  function  of  density,  it 
may  be  possi  le  to  determine  from  the  rotational  temperature  results  if  this 
effect  is  j  sent.  Again,  in  general  the  trend  of  decreasing  Tg  with  increas¬ 
ing  density  seems  counter  to  any  condensation  effect.  A  stronger  point  can  be 
made  once  again  by  referring  to  the  scaling  for  binary  collisions  existing  for 
two  nozzle  diameters  used  at  the  same  (P0d)  value.  Increasing  the  stagnation 
pressure  by  a  factor  of  2-1/2  at  each  (PQd)  value  does  not  destroy  the  scaling. 

It  can  be  noticed  from  the  data  that  for  a  given  diameter,  the  (P0d)  curves 
get  closer  together  as  the  pressure  increases,  and  it  maybe  expected  that  at 
a  high  enough  pressure,  the  trend  may  reverse  itself.  That  is.  at  very  high 
values  of  (PQd) ,  the  curves  would  show  an  increasing  rotational  temperature 
due  to  condensat  .on  effects.  For  a  quick  check  on  this,  a  number  of  data 
points  were  taken  for  (PQd)  =  630  and  930  torr-mm,  corresponding  to  stagnation 
pressures  of  8l5  and  1200  torr  respectively  (with  a  0.75  ran  diameter  nozzle). 

At  x/a  values  greater  than  about  25,  these  data  agreed  approximately  with  the 
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(PQd)  -  480  torr  mm  curve  shown  in  Pig-  16-  At  lower  x/d  values,  the  data  in¬ 
creased  and  rose  above  the  (Pod)  =  ^80  curve.  It  should  be  noted  thet  at  these 
high  pressures,  the  electron  beam  itself  was  scattered,  leading  to  a  very  wide 
fluorescence.  Spatial  resolution  was  destroyed,  especially  with  the  very  small 
nozzle,  and  part  of  the  barrel  shock  may  even  have  been  observed  at  small  x/d 
values  because  the  size  of  the  flow  field  was  very  small  In  that  region-  Thus, 
the  high  pressure  results  are  inconclusive,  except  possibly  at  large  values  of 
x/d. 


Thus,  based  on  the  discussion  of  pressure  dependence,  and  the  fact 
that  the  stagnation  pressures  reported  for  the  present  data  were  less  than  ?.ho 
torr  (i.e.  l/3  atm.).  It  is  felt  that  condensation  plays  no  role  in  the  rotation¬ 
al  temperature  departure  from  equilibrium. 

The  experimental  rotational  temperature  results  can  be  compared  to 
rotational  freezing  values  predicted  by  Knuth  (Kef.  13),  in  which  a  freezing 
criterion  similar  to  that  used  in  chemically  reacting  and  vibrationaliy  relax¬ 
ing  expansions  was  used.  The  freezing  point  for  rotational  energy  Is  defined 
as  the  point  at  which 


= 

Dt 

where  is  the  rotational  energy,  and  Tg  is  the  rotational  relaxation  time. 

If  the  rotational  heat  capacity  can  be  considered  as  constant,  then  Eq.  7  say 
be  rewritten 

■jyp,.  _ 

"Dt:|  "  ^ 

If  the  magnitude  of  DT/l)t  is  either  smaller  or  greater  than  1/tk,  then  the 
rotational  energy  may  be  considered  as  in  equilibrium  with  translation,  or 
frozen.  The  relaxation  time  was  estimated  from  viscosity  calculations 
using  the  Buckingham  ;(6-Exp.)  potential,  and  assuming  a  value  for  ,  the 
number  of  collisions  required  for  rotational  relaxation,  i.e. 

ZR  =  tR  v  (9) 

where  v  is  the  collision  frequency.  Using  a  value  of  Zg  =  5.3  collisions 
from  room  temperature  measurements,  and  employing  the  isentropic  free  jet 
relations,  Knuth  has  predicted  rotational  freezing  conditions  in  free  jets. 

To  compare  the  present  experimental  rotational  temperature  data  with  Knuth 's 
predicted  frozen  values,  a  Tj  freeze  value  is  defined  as  the  asymptote  for  each 
experimental  (P0d)  curve,  i.e. 
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Sxnerirent 


Xnuth 


This  comparison  is  shovn  in  Fig.  30  for  both  the  freeze  point  in  the  free  Jet 
(x/d  freeze),  and  the  rotational  temperature  at  freezing  (Tr  freeze)  ,  It  is 
seen  that  the  rotational  temperature  results  at  low  values  of  (PQd) ,  i.e. 

Pod  <  50  torr-na,  are  in  general  agreement  with  Khuth*  s  predictions  use  a 
collision  number  of  Zr  =  5-3-  However,  at  lower  temperatures,  typified  by  the 
higher  (?Qd)  data,  a  larger  value  of  Zr  is  required  in  the  freezing  criterion 
to  match  the  experimental  results.  Using  a  value  of  Zr  =  10  collisions, 

Smith's  prediction  of  Tr  freeze  agrees  more  closely  with  the  experimental 
values  for  temperature  below  about  20  °K. 

Sufficiently  far  downstream  of  the  orifice  for  a  given  stagnation 
pressure,  the  density  may  become  low  enough  that  translational  freezing  may 
occur,  i.e.  the  free  jet  undergoes  a  transition  to  a  collisionless  free  mole¬ 
cule  flow.  As  used  in  the  present  experiments,  the  electron  beam  probe  will 
not  measure  this  directly,  since  it  is  the  rotational  energy  distribution 
that  is  observed  in  the  beam  fluorescence.  Computations  were  performed  to 
obtain  a  rough  estimate  of  where  translational  freezing  might  occur  for 
some  of  the  present  free  jet  experiments.  Again,  the  procedure  discussed  by 
Xnuth  (Ref.  13)  was  used.  Using  the  predicted  rotational  freeze  point  as 
calculated  by  Khuth,  it  was  assumed  that  the  expansion  proceeded  from  this 
point,  now  with  a  y  =  5/3-  The  "effective"  source  conditions  of  pressure 
and  temperature,  as  well  as  "effective"  source  orifice  diameter,  were  computed 
for  this  expanding  7  =  5/3  gas  using  Khuth' s  relations.  These  equations  were 
obtained  for  a  given  streamtube  by  matching  conditions  at  the  rotational  freeze 
point  for  a  sudden  change  in  7  from  7/5  to  5/3-  With  this  new  expansion  the 
translational  freeze  point  was  predicted,  see  Ref.  13.  In  addition,  the  pre¬ 
diction  for  terminal  (i.e.  freeze)  Mach  number  for  a  monatomic  gas  was  obtained 
from  the  work  of  Anderson,  Andres,  Fenn  and  Maise  (Re?.  47).  A  correlation 
between  terminal  Mach  number  and  ^ource  Knudson  number  was  given  in  Ref.  47 
as  a  linear  function  of  Kh^  _v  -V  vfor  argon  free  jets.  Using  the  "effective" 
source  values  to  obtain  (Kn^) effective >  "the  terminal  Mach  number  was  computed 
for  the  new  expanding  7  =  5/3  gas.  The  results  of  nhese  computations  are 
indicated  below: 
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On  the  basis  of  these  approximate  calculations,  it  can  be  seen  from  referring 
to  Fig.  16,  that  translational  freezing  may  be  occurring  for  the  final  few 
temperature  points  at  each  P0d  (reference  shouldbe  made  to  the  data  obtained 
with  the  5  mm  orifice) .  The  possibility  of  translational  freezing,  then, 
should  be  considered  when  analyzing  the  rotational  temperature  data  in  Fig.  16. 


Based  on  the  experimental  results,  and  the  preceding  dis¬ 
cussions  concerning  the  data,  the  following  conclusions  are  stated: 

1.  Rotational  nonequilibrium  has  been  observed  in  underexpanded  jets 
of  room  temperature  nitrogen.  A  systematic  experimental  survey 
indicates  that  the  rotational  temperature  follows  the  isentropic 
temperature  curve  to  lower  values  as  the  density  is  increased. 

2.  The  electron  beam  fluorescence  technique  may  be  used  at  low  tempera¬ 
tures  to  obtain  rotational  temperatures  in  nitrogen.  The  scatter 
about  the  isentropic  temperature  down  to  about  70°K  is  comparable 
to  room  temperature  scatter,  approximately  -  jflo. 

3.  Effects  of  secondary  electrons  in  contributing  to  a  systematic  error 
in  the  rotational  temperature  measurements  (i.e.  "effective  heating") 
were  not  observed  in  the  present  investigation. 

4.  Possible  effects  of  secondary  electrons  may  contribute  to  the 
apparent  departure  from  a  Boltzmann  distribution  in  rotational  ievels 
measured  during  a  rapid  expansion.  An  experiment  performed  in  a  low 
temperature,  stagnant  gas  should  contribute  more  information  on  this 
phenomenon. 

5.  The  population  in  rotational  levels  becomes  highly  non-Bolt zrnann  in 
the  center  of  shock  waves.  This  effect  becomes  more  pronounced  as 
the  flow  Mach  number  is  increased. 

6.  Measurements  in  free  jets  have  shorn  the  density  to  follow  the  isen¬ 
tropic  curve  in  rotational  nonequilibrium  flow,  following  a  radial 
source  distribution.  The  Mach  disc  at  low  pressures  was  found  to  be 
very  diffuse. 

7.  The  beam  fluorescence  technique  can  be  used  to  obtain  density  pro¬ 
files  through  shock  waves.  The  shock  thickness  data  obtained  in  the 
present  study  agree  quite  well  with  the  results  obtained  in  shock 
tubes. 
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FIG.  1  SCHEMATIC  OF  AN  UNDEREXPANDED  JET  ISSUING  FROM 

A  SONIC  ORIFICE 


001 


FIG.  2  MACH  NUMBER  DISTRIBUTION  ALONG  AXIS  OF  FREE  JET 


TEMPERATURE 


FIG.  3 


TEMPERATURE  DISTRIBUTION  ALONG  AXIS  OF  FREE  JET 


electron  pun  assembly 


FIG.  4  SCHEMATIC  OF  UTIAS  LOW  DENSITY  WIND  TUNNEL 


FIG.  5  SCHEMATIC  OF  VACUUM  SYSTEM  FOR  THE  ELECTRON  GUN 


PIG.  6  ELECTRON  GUN  WIRING  SCHEMATIC 


FIG.  7  SCHEMATIC  OF  SPECTROMETER 


camera  for  flow  visualization 


PIG.  8  TOP  VIEW  SCHEMATIC  OP  EXPERIMENTAL  APPARATUS 


?ection°ofIoEw  dInsi?™  dnnelBLY  mounted  on 
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FIG.  9(c)  PHOTOGRAPH  OF  INSTRUMENTATION  SHOWING  HIGH 
VOLTAGE  SUPPLY  FOR  ELECTRON  GUN,  DEFLECTION 
PLATE  CONTROLS  AND  PHOTOMULTIPLIER  POWER 
SUPPLY. 


FIG.  9(d)  PHOTOGRAPH  OF  TUNNEL  TEST  SECTION  SHOWING  SHOCK 


HOLDER,  FOCUSSING  LENS  AND  5  mm.  NOZZLE  IN  PLACE. 
THE  REMOTE  CONTROLLED  DOOR  AND  BEAM  COLLECTOR 
CUP  ARE  DIRECTLY  ABOVE  AND  BELOW  THE  NOZZLE. 
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FIG.  11 
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FIG.  13 


(continued) 


PIG.  14  R-BRANCH  ROTATIONAL  SPECTRA  OP  N+  0-0  BAND  AT 
VARIOUS  TEMPERATURES  1 


FIG.  14  (continued) 


FIG.  15(b)  LOG  PLOTS  OF  INTENSITY  OF  FREE  JET  SPECTRAL  DATA 
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FIG.  16 


ROTATIONAL  TEMPERATURE  DISTRIBUTION  IN 
UNDEREXPANDED  NITROGEN  JETS 


FIG.  17  ROTATIONAL  TEMPERATC'trHE  MEASUREMENTS  AT 
SMALL  */d. 


RATAL  ROTATIONAL  TEMPERATURE  DISTRIBUTION 
IN  NITROGEN  FREE  JET 


FIG.  19 


PHOTvOMUJLTJPJJIER  DENSITY  ^CALIBRATION  WITH 
ROOM  TEMPERATURE  NITROGEN 


DENSITY  RATIO 


FIG.  20 


DENSITY  DISTRIBUTION  ALONG  CENTERLINE  OF 
NITROGEN  FREE  JETS 
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FIG.  21 


FLOW  VISUALIZATION  PHOTOGRAPHS  OF  SHOCK  HOLDER 
IN  NITROGEN  UNDEREXPANDED  JETS 
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FIG.  22(b)  SHOCK  DENSITY  PROFILE,  Mall 


FIG.  24(c)  LOG  PLOTS  OF  INTENSITY,  M  W  15 


FIG.  25 


COMPARISON  OF  UTIAS  AND  U.  CALIF.  ROTATIONAL 
TEMPERATURE  MEASUREMENTS 


FIG.  26(a) 


TEMPERATURE  DISTRIBUTION  IN  EXPANSION 
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FIQ.  27(a)  SHOCK  WAVE  TEMPERATURE  PROFILE,  M 
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PIG.  27(b)  SHOCK  WAVE  TEMPERATURE  PROFILE,  M 


FIG.  27(c)  SHOCK  WAVE  TEMPERATURE  PROFILE,  M  53  15 


FIG.  28 


VAPOR  PRESSURE  OF  NITROGEN 


TEMPERATURE 


APPROXIMATE  NITROGEN  CONDENSATION  CALCULATION 
FOR  A  FREE  JET 


FIG.  29 
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